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ABSTRACT 

We study the correlations of the shear signal between triplets of sources in the Canada- 
France-Hawaii Lensing Survey (CFHTLenS) to probe cosmological parameters via the mat¬ 
ter bispectrum. In contrast to previous studies, we adopted a non-Gaussian model of the data 
likelihood which is supported by our simulations of the survey. We find that for state-of-the- 
art surveys, similar to CFHTLenS, a Gaussian likelihood analysis is a reasonable approxi¬ 
mation, albeit small differences in the parameter constraints are already visible. For future 
surveys we expect that a Gaussian model becomes inaccurate. Our algorithm for a refined 
non-Gaussian analysis and data compression is then of great utility especially because it is not 
much more elaborate if simulated data are available. Applying this algorithm to the third-order 
correlations of shear alone in a blind analysis, we find a good agreement with the standard 
cosmological model: Eg = <7g(fl m /0.27) 0 ' 64 = 0.79*o ii f° r a flat ACDM cosmology with 
h = 0.7 ± 0.04 (68% credible interval). Nevertheless our models provide only moderately 
good fits as indicated by x 2 /dof = 2.9, including a 20% r.m.s. uncertainty in the predicted 
signal amplitude. The models cannot explain a signal drop on scales around 15 arcmin, which 
may be caused by systematics. It is unclear whether the discrepancy can be fully explained 
by residual PSF systematics of which we find evidence at least on scales of a few arcmin. 
Therefore we need a better understanding of higher-order correlations of cosmic shear and 
their systematics to confidently apply them as cosmological probes. 

Key words: gravitational lensing: weak - cosmology: observations - dark matter - methods: 
statistical 


1 INTRODUCTION 

The statistics of the distribution of matter on large cosmological 
scales, when combined with other cosmological probes, is a pow¬ 
erful tool_todiscriminate_between different cosmological models 
(e.g. lDodelsonll2003l : Laureii s et alfeoilh . Gravitational lensing is 
a technique to assess the mass distribution in the Universe in a way 
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that is independ ent o f the exact nature of dark matter and its dy¬ 
namical state iBartelmann & Schneideij|200ll . for an extensive re¬ 
view). One of the consequences of gravitational lensing is cosmic 
shear, which we stat istically infer from correlations be tween shapes 
of distant galaxies JSchneiderll2006l : [KilbingeJl20l4 for a recent 
review on weak gravitational lensing). The correlations between 
shapes of galaxy pairs give a measurement of the projected matter 
density power spectrum, which in turn constrains the geometry of 
the Universe and the growth of structure. 

The most recent cosmological constraints from cosmic shear 
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repo rted by Fu e t al.l 1 20J4J), F14 hereafter. iKi lbinger et alj 


l 20J_3|f_ Kjtching_etal. J2014 ~). iBeniamin et all d2013h . and 
iHevmans et al .1 J20 1 3l) where the authors analyse the latest data re¬ 
lease by the Canada France Flawaii Telescope Lensing Survey team 
(CFHTLenffl). This lensing catalogue builds upon the Canada 
France Flawaii Telescope Legacy Survey (CFF1TLS), which repre- 


state-of-the-art of gravitational lensing surveys from the ground. 
A preliminary we ak lensing analysis o f the CFFITLS has bee n 


presented earlie r in Hoekstra_etal 


IBeniamin et alj d2007t) ; Fu et al.l 1 20081) . Since then, however, the 


1 20061): Semboloni et all j2006l) : 


CFFITLenS data has significantly improved in terms of the charac¬ 
terisation of the residual systematics and the estimation of galaxy 


reality (Hildebrandt et al. 20121; Erben et al. 20131; Hevmans et al. 

2012 

; van Waerbeke et al. 20131; 

Gillis et al. 20131; Simon et al. 

2013 

;ISimpson et al. 2013;IVelander et al. 20141). 


While current studies mainly focus on the two-point 
correlations in the cosmic shear field, higher-order statis¬ 
tics contain more information, and this can improve con- 


straints on cosmological models (Bernardeau et al. 

19971; 

van Waerbeke et al. 1999; Bernardeau et al. 20031; Takada & Jain 

20031: Schneider & Lombardi 2003|; 

Kilbinger & Schneidei 

2005; 

Schneider et al. 2005!; Berge et al. 

20101: Vafaei et al.l 

201C; 

Kayo et al. 2013:). For the ongoing wide field surveys, such as the 


Kilo Degree Survey (KiDSjfl), the Dark Energy Survey (DEsQ), 
the Flyper Suprime-Cam survey (HSC@), and future surveys such 
as the Large Synoptic Survey Telescope (LSS10) and EuclicQ 
the statistical power of the three-point shear statistics alo ne is 
comparable to th at of two-point shear statistics dVafaei et al.ll2Qlol ; 
Kayo et aHl2013l) . Therefore, prospects on obtaining cosmological 
information from third-order shear statistics are high. 

On the observational side, early attempts to measu re the three- 
point shear statis tics were carried out by IPen et al.l J2003I) and 
I Jarvis et all d2004l) . These two studies were performed using small 
data sets and shape measurement algorithms that are not as ro¬ 
bust as algorithms today. While in both studies a signal was de¬ 
tected, the results were strongly affected by residual point spread 
function (PSF) systematics. More recentlv. lSemboloni et alj <20111 ) 
used high-quality space-based data, the HST/COSMOS dataset, to 
perform a measurement of three-point shear statistics that did not 
show evidence of residual systematics; however the analysis was 
limited to 1.6 deg 2 of the COSMOS data. The latest successful 
me asurements of t hird- order shear statistics have b een performed 
by iFu et al. d2014l) and Ivan Waerbeke et all d2013l) based on the 
CFHTLenS data set. These two different approaches, shear correla¬ 
tion functions and moments in the reconstructed lensing mass map, 
are complementary and are sensitive to different residual systemat¬ 
ics. 

The interpretation of these statistics is still plagued and pos¬ 
sibly limited by theoretical uncertainties. A correct interpretation 
of this signal can only be performed by accurately modelling the 
evolution of the matter bispect rum in the non-linea r regi me. An¬ 
alytical fitting formulae such as IScoccimarro & Couchmanl 1 2001) 
are only accurate at the 10-20% level l van Waerbeke et al J 1200 1 


are only accurate at the 10-20% level Ivan Waerbe ke et all 1 2001 

Semboloni et a l.l201 li fHarnois-Dera ps et alj2012l) . Alternative ap 


proaches based on the halo m o del also have a li mited accuracy 
(IValageas & Nishimichil201 lal lbl; iKavo et al]|2013l) . 

Moreover, other phenomena are expected to affect the mea¬ 
sured signal, such as baryonic physics in the non-linear regime, 
intrinsic alignments of source galaxies and source-lens cluster¬ 
ing. These ph enomena have not ye t been extensively studied and 
are uncertain l Hamana et al] 2002|; Semboloni et al] I2008L120121 : 


iHamois-Deraps et al.ll2014: Fu et al.l 20141). 

Improvements of both the theoretical cosmological models 
and the reduction of systematics in observational data are only 
two pillars of a successful exploitation of the plentiful cosmolog¬ 
ical information in the higher-order shear statistics. The success 
will also depend on realistic models of statistical uncertainties in 
the shear estimators. For this, a Gaussian likelihood is typically 
used in the statistical analysis, such as in F14, whereas at least 
for second-order cosmic s hear statistics there is evidence in favour 
of more complex models | Hartlap eTal. 2009; Keitel & Schneideil 


1201 ll : ISato et al]|201 ll : IWilking & Schneideil 20131) . For this paper, 

we hypothesise that a Gaussian model for the data likelihood of 
third-order shear correlations possibly yields biased results for cos¬ 
mological parameters. We motivate this hypothesis by our obser¬ 
vation in Sect. [4] that the distribution of the estimates in simula¬ 
tions of the CFHTLenS data exhibits a non-Gaussian distribution 
on angular scales of around 10-30 arcmin, violating the assumption 
of Gaussian noise. To test our hypothesis for CFHTLenS data we 
compare the cosmological constraints obtained from measurements 
of the third-moment of the apert ure mass when based on Gau ssian 
versus non-Gaussian likelihoods dSchneider et al J 19981 . 1200^ ). The 
CFHTLenS data is briefly summarised in Sect. [2] Our first analysis 
uses a commonly used Gaussian likelihood as in FI4, whereas the 
second analysis uses a non-Gaussian model. Our estimator of the 
third-order shear statistics is detailed in Sect.[3] For the cosmology, 
we assume a flat ACDM model with the matter density parame¬ 
ter and the amplitude of fluctuations in the matter density field 
ag as free parameters; a flat ACDM model is strongly supported 
by recent constraint s from the cosmic microwave b ackground 
dKomatsu et alj|2009l ; IPlanck Collaboration et alj|2014l) . Based on 
our new technique in Sect. [5] we construct the non-Gaussian likeli¬ 
hood from a set of simulated measurements. 

We present the results in Sect.[6]and discuss them inject. [7] In 
comparison to the two-point systematics analysis of iHevmans et alj 
d2012B . H12 hereafter, we perform new tests for third-order shear 
systematics of CFHTLenS that we present in Sect. l4.2l and in Ap¬ 
pendix 0 


2 DATA 
2.1 CFHTLenS 

The CFHTLS-Wide survey area is divided into four independent 
fields (Wl, W2, W3, W4), with a total area of 154 deg 2 , observed 
in the five optical bands u*, g 1 , r', i', z'. Each field is a mosaic of 
several MEGACAM fields, called pointi ngs. M ore details about the 
data set itself are given in lErben et aD d201 3h . The procedure jor 
the sh a pe measurements u si ng lensfii can be foun d in lMiller et alj 
d2013l) . lMiller et all d2007l) . iKitchjng^taljj 200j|). and the photo¬ 
metric redshifts are described in Hildebrandt et al.1 d2012l) . 


1 http://www.cfhtlens.org 

2 http://www.rcslens.org 

3 http://kids.strw.leidenuniv.nl 

4 http://www.darkenergysurvey.org 


5 http://www.subarutelescope.org/Projects/HSC 

6 http://www.lsst.org 

7 http://sci.esa.int/euclid 
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CFHTLenS: non-Gaussian analysis of 3pt shear correlations 3 


A description of the CFHTLenS shear catalogue, and the residual 
systematics based on the shear two-points correlation function, is 
given in H12. 

For the measurement of the three-point shear statistics pre¬ 
sented in this paper, we use galaxies from the 129 pointings se¬ 
lected by H12, with 0.2 < 2 p hot < 1.3 and i! < 24.7. The mosaic 
for each field has been constructed by merging the single pointings, 
so that overlaps are eliminated, and each galaxy appears only once. 
Each field is projected on the tangential plane cen tred in the middle 
usin g a gnomonic projection (see for example lCalabretta & Greis~eiil 
120001) . Three-point shear correlation functions are measured for 
each of the four fields, i.e., not for individual pointings, using for 
each galaxy the final Cartesian coordinates (x , y) (flat sky approx¬ 
imation), the ellipticity (ei, £ 2 ) and weights w provided by lensfit. 

In order to interpret the shear signal we need to know the red- 
shift distribution of the sources. The redshift probability distribu¬ 
tion function (PDF) of each galaxy in the CFHTLenS catalogues 
is sampled in 70 redshift bins of width 0.05 between 0 and 3. We 
obtain the source redshift distribution of the full dataset by stack¬ 
ing the distributions of all galaxies. However, since each galaxy 
in our sample is weighted according to w when we compute the 
shear signal, we need to weight the PDF of each galaxy to ob¬ 
tain the effective redshift distribution of the sources, p z ( z). This 
technique is explained and tested in lBeniamin et alj 120131) . The fi¬ 
nal redshift distribution, shown in Figure |T] has a mean redshift of 
Zphot = 0.74, and it is sampled in 30 steps between redshift zero 
and 3. 

2.2 Clone simulations of CFHTLenS 

The CFHTLenS clone is a mock survey in which the lensing 
signal obtained from N-body simulations is known, and the 
observational properties, such as galaxy position, ellipticity, 
magnitude, weight, are included such that the clone’s are 
consistent with the data. In this paper, we use the CFHTLenS 
clone for various purposes. A thorough description of the 
dark-m atter-only simulations can be found in lHamois-Deraps et alj 
d2012l) . These simulations have been constructed using the 

WMAP5+SN+BAO cosmology: {f2 m , Ha, , n, h, os} _ = 

{0.279, 0.721, 0.046, 0.96, 0.701, 0.817} iKomatsu et al.l 120091) . 
Starting at an initial redshift of 200, the mass density in the 
simulation box is sampled at 26 different redshifts between 
2 = 3 and zero. The density fields are collapsed along one of 
the Cartesian axes, and the resulting series of planes are used to 
generate shear, convergence, and mass maps inside the light cone. 
The 184 independent line-of-sights cover an area of 12.84 deg 2 
each and have been populated with sources using the same redshift 
distribution and galaxy density as in the CFHTLenS observations. 

The clone does not include density fluctuations larger than 
the simulation box, and this is known to affect the covariance es¬ 
timated from these simulated maps. These missing supe r-survey 
mode s propagate in many ways in the shear covariance jLi et alj 
120141) . Our measurement, however, is very weakly impacted by 
these for two reasons. First, our shear aperture statistics uses only 
aperture scales up to 30 arcmin, which is well below the max¬ 
imum usable scale of 70 arc min. T his makes the Finite Sup¬ 
port effect described in lHarnois-Deraps & van Waerbekel d2014l) 
at most a 5 percent deficit on the error bar of the largest an¬ 
gles. Second, we must also examine the contribution from the 
Halo Sampling Variance, Beat Coupling and Dilation, which 
causes the small scal e clustering variance to b e under-estimated 
dHamilton et al.l I 2 OO 6 I : iRimes & Hamiltonll20 06). As summarised 



0 12 3 4 
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Figure 1. Source redshift distribution p z (z) obtained by stacking the 
probability densities of all galaxies with 0.2 < z p hot < 1.3, i' < 24.7, 
weighted by the lensfit weight w. This distribution is used here to construct 
a forecast of the cosmological shear signal for a WMAP5 cosmology. 

in lHarnois-Deraps & van Waerbekel (120141) . simulation boxes of 
500 /i _1 Mpc miss about 90 percent of the non-Gaussian part of the 
variance at 2 = 0. We can expect that the clone misses even more 
due to the smaller size of the box. However, weak lensing projects 
many scales onto the same angular measurement, which dramati¬ 
cally decreases the non-Gaussian contribution to the error bar. In 
this case, the most important contribution to the super-sample co- 
variance comes from the Halo Sampling Variance , which peaks at 
small scales. As argued in Kilbinger et al.l ( I2013I) . this causes the 
small scale covariance to be under-estimated by less than 10 per¬ 
cent, hence we do not explicitly correct for super-survey modes in 
the clone. 

3 THREE-POINT SHEAR STATISTICS 
3.1 Definition of statistics 

In this section, we briefly outline the relation between cosmic shear 
and the statistics of the matter density fluctuation. The matter den¬ 
sity field at comoving position x an d at redshift 2 is p(x; z ) = 
p{z)[ 1 + <5(x; 2 )], where 8 is the density contrast and p(z) the av¬ 
erage density at redshift 2 . The power spectrum P(k; 2 ) and the 
bispectrum B{k\, k 2 , kz\ 2 ) are defined by the correlators 

^5(fci; z)5(k 2 \ z)^ = (27r) 3 (5o(fei + k 2 )P{kr,z) (1) 

and 

(ft{k\-, z)S(k 2 ; z)<5(fc 3 ; 2 )^ (2) 

= (27r) 3 fo(fci + k 2 + k 3 )B(k 1 ,k 2 , k 3 ; 2 ) , 

of the Fourier transform 8(ki ; 2 ) of the density contrast <5 for the 
3D wave number ki at redshift 2 . 

We apply a flat-sky approximation in the following. Let x_l 
be the two dimensional vector that we obtain by projecting x onto 
the tangential plane on the celestial sphere that is defined by the 
line-of-sight direction; the x- and y-coordinates of x± are *, y re¬ 
spectively. The complex shear y = 71 + ry 2 and the convergence 
k at angle s are both functions of the second-order derivatives of 
the gravitational potential 4>(x±', x) at X± := s /k(x) with /k(x) 
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being the comoving angular diam eter distance at radial distance 
X = Ixl- Following ISchneider et alj d 19981) we find at an; 
the tangential plane: 


by: 


k(s) 

7i( s ) 

72 (s) 




3Ontf 0 2 f x \ ,g(x)Mx) 


2 c 2 


f 


d X J ""-p id'dytfa X ) 


; in 

r 0 (s,t') 

= (7( r )7( r + s )7 (r + t )e _6 ““) , 

(11) 

(3) 

TiM') 

= <^7*( r )7( r + s )7( r + t )e _2 ‘“) 1 

(12) 


r 2 (s,t') 

= ^7( r )7*( r + s )7( r + t)e _2 *“) , 

(13) 

(4) 

r 3 (*,t') 

= ^7( r )7( r + s )7*( r + t)e _2 *“) , 

(14) 


where s = 

V ss* of s and t/ = ts*/s. In equations ill lb d 14b. the 


(5) 


where Ho is the Hubble constant; fi m is the matter density param¬ 
eter; c is the vacuum speed of light; a is the cosmological scale 
factor; \h is the size of the Hubble horizon; andp x (x)dx in 


ensemble averages are performed over triangles invariant under 
translation and rotation. These triangles can hence be characterised 
by three real-valued variables s and t' = t/ + it 2 . The component 
To is invarian t under a parity transformation, whereas 17 , 2,3 are 
not dSchneider & Lombard ill 20(13!) . 

For our analysis, we integrate the natural components to obtain 
an alternative third-orde r statistic o f shear, the third moment of the 
aperture mass ^Schneider et all 1998h . 


,w : = fWfe/ <«> 

describes the distribution of sources per distance interval dx, or per 
redshift interval in the case of p z ( z)dz. 

In the weak lensing regime of cosmic shear, we find 
1 and |'y| 1 such that the observable galaxy ellipticity 

e = ei 4 - ic 2 becomes an unbiased estimator of the galaxy shear 
to first approximation, i.e., (e) = 7 . In addition, by assuming that 
galaxies are randomly oriented intrinsically, the Fourier transform 
of the angular correlation (ejej) = ( 7 , 7 *) = (KiKj) between the 
ellipticities of pairs of galaxies i,j can be interpreted as a direct 
measure of the matter power spectrum in projection on the sky, 


P ri (£) = P 1 (£) 


9^Pq 4 

4c 4 



dx 


ff 2 (x) 

° 2 (x) 


P 



By 7 * we denote the complex conjugate of 7 . The angular power 
spectrum P K (£) of k is defined by 


(k(£i)k(£ a)> = (2tt ) 2 MTr + ^)P.{ 11 ) (8) 


for the angular wave number £. 

A similar relation exists between the angular bispectrum B K 

of K, 

(k(£i)k(£ 2 )k(£ 3 )) = ( 2 tt) 2 Sd(£i+£ 2 +£ 3 )B k (£i, £2,£3) , (9) 


attainable through the correlation of three galaxy ellipticities, and 
the projected matter bispectrum 


B k (£i,£2,£3) 


27-HpfIm f g 3 (x) 

8 c 6 J X /k(x)® 3 (x) 

V/k(x) , /k(x)’/k(x)’ W ) 


( 10 ) 


jBemardeau et alj 1 19971 ; ISchneider et alj 1 19981) . In contrast to 
second-order statistics, four independent correlation functions of 
the third-order statistics can be defined. We utilise the representa¬ 
tion of th e correlator in te rms o f its n atural components as advo¬ 
cated by I Schneider & Lombardi 0003 ). Let r, r + s, and r +1 be 
the angular positions of three galaxies, while a denotes the angle 
between the x-axis and s. The natural components are then given 


M ap (0) = jd 2 rUe(\r\)K(r) , (15) 

for an aperture filter Ue(r) and a triplet ( 81 , 62 , 83 ) of aperture 
radii, 


< 114 ,( 0 !, 02,03)) 



d 2 r,d 2 /,r„,(r ; )LTL 

(27t) 2 


(16) 


(k(£i)k(£ 2 )k(£3)) 


( Pen et al . 2003l : ljarvis et all2004l : lschneider et al.l2005l) . Roughly 
speaking, M ap ( 0 ) is the Ue -smoothed convergence field at r = 0 , 
and <Mf p ( 0 i, d 2 , $ 3 )) is the [/g-smoothed version of the correla¬ 
tor (k(£i)k(£ 2 )k(£ 3 )) , which is a measure of the projected matter 
bispectrum in Eq. m In particular, by choosing an exponential 
aperture filter as in Ivan Waerbekel J19980 . 


Ue(r) 



r \ r,2 

_ ] p 2 t e > 

20 2 / 


07) 


the relation between three-point correlation functions and the third 
moment of (Aff p ) is relatively simple. Note that our filter has been 
rescaled by 2y/28 in comparison to van Waerbekel J 1 9981) . Our filter 
peaks at £ = 1 / 2/8 « 4862 (0/l ') _1 in wave number space. 


3.2 Model of the matter bispectrum 

The equations & Gil, and GD establish an explicit relation be¬ 
tween B(k\, k, 2 , fe; z) and the third moment <Mf p ). Thus, in or¬ 
der to predict <M 3 p ) we need to model the bispectrum of mat¬ 
ter fluctuations. It is important that the bispectrum model cap¬ 
tures the mode coupling beyond perturbation theory, since our 
CFHTLenS measurements are probing the non-linear regime of 
k ~ 0 .1 — 10 /tMpc^ 1 . To this end, we employ the bispectrum fit 
of lScoccimarro & Couchmar] 1200ll SC01 hereafter). The work of 
SC01 produced an analytical fit to the non-linear evolution of the 
bispectrum based on a suite of cold dark matter N-body simulations 
that was available at that time. The accuracy of this fit is limited by 
the accuracy of the N-body simulations in their study and by the 
fact that the effect of baryons on the small-scale clustering of mat¬ 
ter is not included. In short, the fit of SC01 consists of a refinement 
lowest-order perturbation theory, 

B(ki, k 2 , k 3 ; z) (18) 

= 2F(ki,k 2 ; z)P(kv, z)P(k 2 \ z) + cycl. 
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where cycl. indicates a cyclic permutation of the indices. SC01 
express the non-linear extension of the mode coupling factors 

F(ki, k, 2 \ z ) as 


F(ki,k 2 ,z) ~ -a(n,ki', z)a(n,k 2 ; z) 


(19) 


+ 


1 fci • &2 /fci , k 2 

2 kik 2 
2 / fci • k2 \ 2 


(.Tf + ki ; z)b(n, k 2 ; z) 


+ y( l ik2 2 ) c(n,kv,z)c(n,k 2 ;z) 


where the coefficients 

1 + a 8 ( a )- 0 - 2 [0.7Q3(n)] 1/2 (g[z]/4) Tt+3 - 5 


a(n, fc; z) = 
b(n,k\z) = 
c{n,k~,z) = 


1 + (g[z]/4)"+ 3 - 5 
1 + 0.4(n + 0.3)q[z]" +3 
1 + q[z] n + 3 - 5 

1 + 4.5/[1.5 + (n + 3) 4 ](2q[z]) n + 3 


1 + ( 2 q[z ]) n + 3 - 5 

have been fitted to the N-body simulations with 

4-2" 


Q 3 (n) = 


1 + 2 "+! ' 


(23) 


Here as(z) is the standard deviation of matter density fluctua¬ 
tions within a sphere of radius 8/t^ 1 Mpc linearly devolved from 
zero to redshift z, and n is the spectral index of the primordial 
power spectrum. The time dependence of F(ki, k 2 ; z) is given by 
both the evolution of erg(z) and the function q[z\ = fc/fc n i[z] where 
47 rfc n i[z] 3 Pij n (fc n i[z]; z) = ldefinesthe wave number k n \ [a] ofthe 
non-linear regime at redshift z; Pn n (fc; z) denotes the linear matter 
power spectrum. 

SC01 showed that this app roximation is accur a te to within 
15% up to k of a few foMpc -1 . Ivan Waerbeke et aD d200lh com¬ 
pared the third-order moments of the projected density field mea¬ 
sured directly on simulated k maps with predictions obtained using 
the fitting formula. They f ound a sim ilar ac curacy . In agreement 
with these previous results. [Semboloni et al.l J201 If) found that this 
approximation systematically underestimates (M ap ) on small an¬ 
gular scales. A different approach to com pute the bis pectrum has 
been recently suggested by IValageas & Nishimichil j201 lal fbh. It 
uses a combination of perturbation theory and the halo model. This 
approach is promising but its performance depends on the accu¬ 
racy of the halo-model which is in general still limited. Moreover, 
none of these approximations accounts f or the poten tially large ef¬ 
fects from baryonic physics I lSemboloni et alj|2012l ). Overall, the 
accuracy of bispectrum predictions is therefore still an open issue. 
Current models cannot claim an accuracy better than ~ 20% which 
we include in our cosmological analysis. Consequently further im¬ 
provements in modeling are necessary in the future, which is be¬ 
yond the scope of this work. 

We no te here that the c oeffici ents of SC01 have recently been 
updated by iGil-Marin et all J2012h . Our analysis does not include 
this update as the corrections are smaller than the 20% model er¬ 
ror that we include in our analysis, and as such this update would 
not impact our results. Moreover, as non-li near power spectrum 
P(k\ z) in Eq. ( 1 1 8t we use the model of Smith et ail(l2003l) with the 
tr ansfer function o f lEisens tein & Hul (1998). As recently reported 
in lHarnois-Deraps et al 1 20141) . this model lacks power on small 
scales in comparison to N-body simulations. As shown in Sect. 14. 1 1 
however, this bias is negligible on the angular scales that we exploit 
for our analysis. 


3.3 Estimators of the natural components 

The measurement of the third-order moment of the aperture mass 
statistics, (M 3 „ ( 61 , 62 , 63 )), is pe rformed using the same proce¬ 
dure described in ISemboloni et all J201 ll) and in the recent F14 
analysis (see their Sect. 2.3). This procedure consists of recon¬ 
structing (M 3 p) by numerical integration of estimates of the natu¬ 
ral components T,. 

For the estimators of the natural components equations (DU¬ 
GS. we bin the products of three source ellipticities of similar 
triangle configurations. One possible choice is to bin the trian- 
les according to their values of (s,t'i,t' 2 ). However. Ijarvis et alj 
20041) suggested a more suitable binning scheme. Given three sides 


(20) 

5, 

t, di with s < t < di, the following variables are defined: 


(21) 

d 

— 5 , ffmin <1 d "C dmax , 

(24) 

u 

= s/t , 0 < U < 1 , 

(25) 


V 

= , -1 <V < 1 . 

(26) 

(22) 




Assigning a sign to v keeps track of the triangle orientation: if 
t' x = t x (that is s x > 0 , i.e., the triangle is clockwise oriented) then 
v > 0, otherwise v < 0. The limits d m i n and d max are the mini¬ 
mal and maximal lengths for the smallest triangle side s and define 
the range on which the natural components are sampled. For the 
CFHTLenS fields, d m i n is set to 9 arcsec in order to avoid bias from 
close galaxies pairs with overlapping isophotes. The maximum sep¬ 
aration is set to d max = 400 aremin which means that (M ap ) can 
be reconstructed up to angular scales < 100 aremin. We compute 
the correlation function of the galax y triplet s by_a tree-code ap¬ 
proach, similar to the one suggested bv IZhang & Pen! d2005f) . In our 
implementation of the tree code, we require (Si + S 2 )/D < 0.1 as 
criterion for stopping a deeper search into the tree; Si and S 2 are 
the sizes of any pair of tree nodes belonging to a node triplet, and 
D is the distance between the centres of the nodes. 

In the case of the CFHTLenS catalogue, it is necessary to ac¬ 
count for a multiplic ative correction fac tor m(izsN, r ga i) assigned 
to each galaxy (H12. lMiller et al.ll2013h . The correction factor de- 
pends on the galax y signal-to-noise i/sn and size r ga i. According to 
Mill er et al j d2013l') . for an average shear ( 7 ,) the corrected estimate 
of the i-component is given by 


(7i)cal = 


(7i) 


, with * = 1,2 


(27) 


/cal (1 + m) 

where (...) indicates ensemble weighted averages forn ga i galaxies 
defined using the lensfA weights w, namely 

-1, 


( n gal 

Jfuh 

a=l 


n ga,l 


.w a 


(28) 


/"gal \ - 1 "gal 

(1 + m ) = EH a(l + W a (//SN, / gal )) • (29) 

\a=l / a=l 

The extension of this calibration scheme to the three-point shear 
statistics is straightforward: 

) ] CijaCybe/CjctnatUbmc 

(7*7j7fc)cai ^ u> a w b w c (l + m a )(l + mb)(l + m c ) 

abc 

where the sum is over all triplets a, b, c belonging to a bin and 
i,j, k,= 1, 2 denote the shear components. This correction is ap¬ 
plied to our practical estimators of the natural components, 


r 0 (s,O = 


E 


a,b,c 


w a WbW c e a eb€ c e 


-6io- b 


Erin WaWbWc(l + m a )( 1 + m b )(l + m c ) 


(31) 
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and 

, Er© c w a w b w c e;e b ece -2l “ b 

ri(s,t') = =- a ’ b ’ c , ------ . 

ErET w a w b w c (l + m a )( 1 + m b )(l + m c ) 

(32) 

We obtain the remaining two other estimat ors r 2 . 3 (s,t / ) b y cyc lic 
perm utations of the triangle parameters l lSchneider & Lombardil 

3.4 Estimators of modes of the aperture statistics 

In practice, we encounter four distinct modes of the aperture statis¬ 
tics due to the possible presence of B-modes in the shear field. 
Known sources of B-modes are systematics due to residuals in the 
PSF correction or intrinsic alignments b etween the sources (e.g., 
iHeavens et alfeOO&lKitching et al .1120 1 2t) . Even without these sys¬ 
tematics higher-order effects such as lens coupling and corrections 
beyond the Born approximation give rise to B-modes, however 
much smaller than the amplitude of the E-mode s dCoorav & Hul 
l2002l:lHirata & Seliakll2003l:ISchneider et alJll998h . 

For the aperture mass of an aperture centred on the origin r = 
0, the formal separation between E- and B-modes is given by 

M{ 8 ) = M ap (0) + iM ± ( 8 ) (33) 

= - J d 2 rQg(\r\)~i{r)e~ 2lc1, , 

where the imaginary part M±( 8 ) is the B-mode of the aperture 
mass, the real part M ap (0) is the E-mode of Eq. d 15b . <j> is the polar 
angle of r, and 

Qg(x) = ^^2 J dssUeis)^ - Ug{x) (34) 

is the filter of the shear field that corresponds to Ug. We denote by 
M*( 8 ) the complex conjugate of M( 6 ). To lowest order in Scp/c 2 , 
gravitational lensing can only produce E-modes, hence M± = 0 in 
this regime. In this case, the third-order moment (M 3 p (0i, 62 , 83 )) 
in Eq. ( fT6t is given by the average (M(0i)M(0 2 )M(0 3 )} over all 
available aperture positions. Generally, however, we have E/B mix¬ 
ing inside the correlator, giving rise to different modes of the statis¬ 
tics: the EEE, EEB, EBB, and BBB mode, 

EEE : (M a 3 p (0!,f9 2 ,0 3 )) (35) 

= i Re «M 3 (0 1 ,02,0 3 )) + <M 2 M*(0 1 , 02,0 3 )) 

+ (M 2 M*(0 3 ,0 1 ,0 2 )) + <M 2 M*(0 2 ,0 3 ,0i))) 

EEB: <M a 2 p M x (0 1 ,02,0 3 )) (36) 

= i Im ((M 3 (0 1 , 0 2 , 0 3 )) + (M 2 M*(0 3 , 0i, 0 2 )) 

+ (M 2 M*(0 2 , 0 3 ,0i)> - ( M 2 M * (01,02, 0 3 ))) , 

EBB: <M ap Mi(01,02,0s)) (37) 

= I R e(-(M 3 ( 0 i, 02 , 0 3 )) + <M 2 M*(0i,0 2 ,0 3 )) 

+ (A/ 2 M* (0 3 ,0l,02)) -(M 2 M*(0 2 ,03,0i))) , 

BBB: (Ml (01,02,0s)) (38) 

= i Im ((M 3 (0 1 , 0 2 , 0 3 )) + (M 2 M\ 8 1 ,02,0 3 )) 

+ <M 2 M*(0 3 ,0i,02)) + <M 2 M*(02,03,0i))) ■ 



1 10 

d(arcmin) 


Figure 2. Measurement of the equilateral for the noise-free data of 

the CFHTLenS clone. The signal has been divided in EEE (black solid 
squares), absolute values of the two parity modes BBB (blue) and EEB 
(red), and the B-mode EBB (green). Error-bars represent the error on the 
average computed using the full suite of clone mo ck catalogues. For com- 
parison, we show WMAP5 predictions done using iSmith et al.l 12003 ) and 
the fitting formula SC01 in Sect. 13.11 

Note the order of the arguments (0i, 82 , 83 ) in the previous equa¬ 
tions. The actual cosmological signal EEE is the focus of this study. 
The modes EEB and BBB, or parity modes, are only present for 
parity violation (jSchneiderll2003 D. whereas the B-mode EBB may 
be the product of intrinsic alignments of the galaxy ellipticities or 
an indicator of PSF systematics. To obtain the different modes, we 
estimate the two statistics (M 3 ) and (M 2 M *) from the measured 
natural components as in F14 (their equations 16 and 17). 

Ideally, all four modes can be separated perfectly if the natural 
components are measured with infinite resolution and over the en¬ 
tire sky. However, the observed shear fields are finite, incomplete, 
and they are only sampled at the positions of sources. For second- 
order statistics, specially designed filte rs can be used to do this job; 
see for example I Schneider et al J OOlOf) and references therein. For 
the third-order statistics in this study, on the other hand, this sep¬ 
aration is currently not perfect. However, all these effects can be 
quantified and are much smaller than the expected signal. In addi¬ 
tion, given the large errors in our measurements these inaccuracies 
are not important for this study. Therefore separating the signal into 
E-. B-, and parity components is nevertheless a very effective way 
to assess the cosmological origin of measured shear statistics. 

4 SIMULATED DATA VERSUS MEASUREMENTS 
4.1 Clone Simulation 

We measure the three-point shear signal on the 184 simulated 
12.84 deg 2 lines-of-sights from the CFHTLenS clone. For this 
purpose, we perform two separate measurements. Both sets of mea¬ 
surements are used for different purposes. First, for the noise-free 
sample, we analyse shear catalogues that assume intrinsically round 
sources. Second, for the noisy sample, we add ellipticity noise with 
an amplitude similar to the CFHTLenS data. The noisy sample is 
attained by adding a Gaussian random value with an average of zero 
and a ID dispersion of cr e = 0.28 to a noise-free sample where this 


© 0000 RAS, MNRAS 000, 000-000 















































CFHTLenS: non-Gaussian analysis of 3pt shear correlations 7 


value is the measured ellipticity dispersion of the CFHTLenS in¬ 
cluding both intrinsic ellipticity dispersion and measurement noise; 
we do not include intrinsic alignments. For each line-of sight, we 
measure the natural components r, using the same binning and the 
same criteria we use on the CFHTLenS data. From the natural com¬ 
ponents, we then compute the third-order moments of the aperture 
mass by a numerical integration. 

Figure [2] verifies our method of computational analysis, sup¬ 
porting previous findings that the SC01 matter bispectrum is a rea¬ 
sonable fit to a dark-matter only simulation of the standard cos¬ 
mological model. Using the noise-free sample the figure displays 
the mean and standard error, due to cosmic variance, of the equi¬ 
lateral (M®p(0i = 62 = #3)) = (Ma P (0)) obtained by averaging 
the 184 clone fields. The signal has been separated into E-, B- 
and two parity modes. As mentioned before, due to E/B mixing of 
insufficiently sampled fields we expect a small level of B-modes in 
the estimators even for the pure cosmic shear fields as the ones in 
the simulation. We find that the BBB and EEB modes are consis¬ 
tent with zero, whereas the EBB mode exhibits a positive signal, 
albeit almost two orders of magnitude lower than the cosmologi¬ 
cal signal EEE. A possible explanation for this (negligible) EBB 
contamination, apart from the mixing, might be a numerical resid¬ 
ual due to our numerical transformation from natural components 
to the aperture statistics. For reference, in the figure we also show 
the prediction of (M® p ) for a WMAP5 cosmology and our SC01 
model of the matter bispectrum. Moreover, we include the redshift 
distribution p z (z) of sources as shown in Sect. 12. II 

The agreement of the prediction and the measurement for the 
EEE signal is good overall. There is a discrepancy of about 10% at 
around 5 arcmin increasing to 30% at 1 arcmin. The discrepancy 
may be explained by the limited resolution of the simulations and 
the limite d accuracy of our model of t he no n-linear matter power 
spectrum dHamois-Derans et alJl2012Ll2014l) . In addition, there is 
a comparable inaccuracy in the non-linear regime due to baryonic 
physics that is not accounted for in the clone. We therefore take 
a conservative stand in our analysis and utilise only measurements 
on scales larger than 5 arcmin. 

In contrast to our rel atively good ag reement between simula¬ 
tion and analytical model. IValaeeas et all d2012h find that the SC01 
model underestimates the power of a simulation by roughly a fac¬ 
tor of two at 5 arcmin of (Mf p ) (their Fig. 2; lower left panel). 
However, the authors in this paper employ a polynomial filter that 
peaks at 0 DO i v ~ 5/1, whereas our Gaussian aperture filter peaks 
at 9 « s/2/t dSchneider et al.lfl998l : ICrittenden et alJl2002h . This 
means that the corresponding scales of the polynomial filter are 
0 po iy >17.6 arcmin for 9 > 5 arcmin, where the agreement of 
SC01 is much better compared to the simulation. All the same, it 
is still possible that our good agreement between simulation and 
model is a coincidence that results from an offset of the clone 
power due to the limited resolution and volume of the simulation. 
In order to further investigate this, we reproduce the Fig. A1 of our 
companion paper F14, see our Fig. [3] In this figure, we compare 
the clone measurements to a suit of recent, more accurate mod¬ 
els of the matter bispectrum discussed in the literature (see figure 
caption). All these models agree with the clone within ~ 10% 
for angular scales larger than or equal 5 arcmin (Gaussian aperture 
filter). In addition, our model SC01+S03 is below the amplitude of 
all recent models by about 10%-20% for scales larger than 5 ar¬ 
cmin (bottom black line). We therefore conclude that (i) analytical 
models of the matter bispectrum on the scales considered here root- 
mean-square vary within ~ 20%, and (ii) the bispectrum power in 
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Figure 3. Data points of the equilateral (Mf p ) in the clone (open circles) 
in comparison to different analytical descriptions of the matter bispectrum 
(lines). The models are combinations of a fitting formula for the bispec¬ 
trum with different prescriptions of the n on-linear matter power spectrum ; 
SC0 1: IScoccimarro & Couchmanl 1200 ll): GM12: Gil-Marin ct all j2012l); 

H13: Heitmann et affl2014l) : T12: lTakahashi et al.|j2012l) : S03: lsmith et alj 

I2OO3I1 . We use SC01+S03 for the scope of this study (bottom black line). 
The figure is a reproduction of Fig. A1 in FI4. 


the clone is sufficiently well described by SC01 on the angular 
scales that we consider by a cosmological analysis. 


4.2 CFHTLenS measurement 

In this section, we present our measurements of (Mf p ) in the 
CFHTLenS data. To start, in the left panel of Fig. [4] we evaluate 
the difference between (Mf p ) as measured using all the 129 point¬ 
ings to the same quantity measured excluding 9 fields that fail our 
new residual systematics test in Appendix [A] Our systematics test 
refines the tests done in H12 for third-order shear data. In the fol¬ 
lowing, we refer to the 120 remaining fields as pass fields. These 
sample is utilised to constrain cosmological parameters. For Fig. [4] 
we calculate the error bars by rescaling the standard error of the 
mean in the 184 noisy clone simulation by the factor VA ra tio, 
where A ra tio = 0.12 is the area ratio between the pass fields and 
the simulation. The amplitude difference between the two measure¬ 
ments of (M® p ) with and without the rejected 9 fields are smaller 
than the statistical error. Nevertheless the two measurements are 
performed essentially from the same sample, therefore the change 
in amplitude at aperture radii >15 arcmin may be indicative of 
residual PSF systematics on these scales. The right panel of Fig. 
[f] shows the CFHTLenS measurement of (Mf p (0)) for the pass 
fields in comparison to the predicted signal for a WMAP5 cosmol¬ 
ogy, and again the measurement from the clone. 

In Fig. [5] we compare our measurement of the EEE mode to 
the other modes of (Mf p (0)), and we possibly find more evidence 
for PSF systematics on scales of a few arcmin and below. For scales 
9 > 3 arcmin, the EEE signal is roughly one order of magnitude 
larger than the B- and parity modes whenever it is not consistent 
with zero. Meanwhile, the amplitude of the B- and parity modes are 
for 9 < 5 arcmin clearly larger than those measured on the idealis¬ 
tic data in the clone. Fig. [2] where typical amplitudes are always 
smaller than ICC 9 . Possible origins of the increase of the system- 
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Figure 4. Left panel : the black solid triangles show the equilateral computed using all the 120 CFHTLenS fields passing the systematic tests by 

H12 and our additional test in the Appendix I aI The red squares show the same signal measured using all 129 fields that passed H12 tests only. The error bars 
have been computed using the clone. Right panel : display of the measured from the 120 pass fields together with the WMAP5 prediction and the 

average signal from the clone. 
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Figure 5. Left panel: comparison of the CFHTLenS EEE mode of (black data points) to the absolute values of EEB, EBB, and BBB modes (red, 

green, and blue lines). The error bars have been computed using the clone. The dashed black line is a WMAP5 prediction. Right panel: same quantities as in 
the left panel after excluding elliptical galaxies, i.e., sources with BPZ Tg ^ 2. 


atics indicators can be intrinsic alignments (IA) and PSF residuals, 
both of which are not accounted for in the clone. In principle, 
IA can generate B-modes, but the correlations in the source ellip- 
ticities should remain parity invariant. The presence of both parity 
modes EEB and BBB signals hence indicates that IA alone can¬ 
not explain the systematics in the data. This hints to an imperfect 
PSF correction in the shape measurement that may be relevant for 
(•Map) - This further supports our decision to reject measurements 
below 5 arcmin in the cosmological analysis. 

Nevertheless IA co uld also affect the measurement of our cos¬ 
mological shear signa l l lMerkel & Schafed I20l4 IValageasI 120141 : 
ISemboloni et alj|2008l) . For this reaso n, we try to quan t ify the im¬ 
pact of IA. Using N-body simulations. lSemboloni et alj ( 120081) find 
that mostly the IA of early type galaxies contaminates three-point 


shear statistics. In their model intrinsic ellipticities of galaxies are 
given by the ellipticity of the parent matter halo. Early type galax¬ 
ies are perfectly aligned to the halo, whereas spiral galaxies have 
a random misalignment to their parent halo. In surveys compara¬ 
ble to CFHTLenS wide, the IA amplitude of the EEE signal due to 
intrinsic-intrinsic-intrinsic (III) correlations becomes strong com¬ 
pared to the cosmic shear signal for aperture radii below a few 
arcmin and grows to a comparab le amplitu de at ~ 1 arcmin. By 
considering tidal torque theory iMerkel & Schafei] d2014l) find that 
third-order correlations between the intrinsic shapes of sources and 
shear (GGI) are negligible with respect to the III correlations for 
our angular scales of interest; the absolute amplitude of III and GGI 
correlations is highly uncertain in their model though. Using the as¬ 
sumption of a linear relation between intrinsic shapes and the local 
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matter density as well as hierarchical clustering, IValageasl i2014l) 
finds that IA contamination is of the level of typically about 10 % 
for different source redshifts which is ubiquitous in our data. This 
figure includes a bias due to lens-source clustering which arises 
because of a correlation of source number density and the matter 
density (treated as linear in the model). In summary, these pieces 
of research imply that IA alignments should make only a small con¬ 
tribution to the total EEE signal for 6 > 5 arcmin but may be sub¬ 
stantial for 9 ~ 1 arcmin. 

Using the conclusions from lSemboloni et alj ( l2008h . it is pos¬ 
sible to evaluate the impact of III in our analysis by measuring the 
three-point signal with and without elliptical galaxies. The effec¬ 
tive removal of ellipticals is achieved by selecting galaxies with 
classification Tb > 2 from the Bayesian Photome tric Redshift Es¬ 
timation (BPZ: lBenitej2000l ; lHevmans et al.l20l 3l) . This decreases 
the effective number of sources by about 20 percent. With this se¬ 
lection the change in the EEE signal is within the error bars, as we 
show in the right panel of Figure[5] The EBB, EEB and BBB com¬ 
ponents qualitatively retain their behaviour compared to the full 
galaxy sample: their amplitude is large at small aperture radii and 
falls off quickly. The amplitude itself, however, changes in a way 
that is difficult to reconcile with IA: the BBB signal increases sig¬ 
nificantly at about 3 arcmin rather than decreasing as anticipated. 
Our result therefore suggests that (i) the systematics are associated 
to residual PSF systematics which affect late-type galaxies more 
than early-type galaxies and (ii) that the precision of the measure¬ 
ment is still not sufficient to set constraints on the IA from the 
three-points statistics alone. For reason (i), we decide not to ap¬ 
ply a morphological cut for our cosmological analysis. Hence our 
final catalogue is composed of 120 pass fields with the original se¬ 
lection cuts previously described in Section[2] The final selection of 
pointings correspond to a total effective area of roughly 100 square 
degrees, with masked regions taken into account. 

Figure [ 6 ] shows, for this final catalogue, the measurement of 
the full 62 , # 3 )) in comparison to the WMAP5 predic¬ 

tions and the clone measurement. As before with the equilat¬ 
eral (Mf p ( 0 )) we find a reasonable agreement also between the 
full statistics of the measurement, WMAP5 prediction, and the 
clone. For 0; < 5 arcmin, the SC01 model falls slightly below 
the clone, on larger scales the CFHTLenS measurement appears 
to be below the prediction; the statistical errors are, however, rel¬ 
atively large in this regime. In the next section, we use this vector 
of data points for 9i ^ 5 arcmin to infer cosmological parameters 
through a likelihood analysis. 


4.3 Evidence for a non-Gaussian likelihood 

The matter density field obeys non-Gaussian statistics for smooth¬ 
ing scales below the typical size of large galaxy clusters, and it 
asymptotically approaches Gaussian statistics towards larger scales 
and higher redshifts for a Gaussian primordial density field. Be¬ 
cause of this non-Gaussian nature on small physical scales we ex¬ 
pect that measurements of the three-point statistics of cosmic shear 
will exhibit a distinct non-Gaussian distribution on small angular 
scales, if these scales are not dominated by the shape noise of the 
sources. In order to investigate whether the estimator of (M® p ) 
shows any signs of non-Gaussianity, we study its distribution in 
184 realisations of the clone for both the noise-free and the noisy 
samples. For this purpose, we compute estimates of (Mf p ) for the 
same combination of aperture radii as for the CFHTLenS data. The 
resulting 184 simulated data vectors represent the likelihood of ob¬ 
taining a value of on a 12.84deg 2 survey (correspond- 



<Mj(0)> <mJ(9)> 


Figure 7. Top panels', frequency of values of the equilateral (ATf p (0)) from 
184 lines-of-sight from noisy version of the clone (12.84 deg 2 ). In the 
left panels, we show the result for 9 = 3.438 arcmin, while in the right 
panels we show the distribution for 6 = 30 arcmin. Black dashed lines 
show the average values, while the solid red lines, indicate the best-fitting 
Gaussian. The black solid line indicates the zero value. Bottom panel', the 
same as the top panels for the noise-free version of the clone. 
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Figure 8. Top panel', median signal and the ±25% quartiles obtained from 
the 184 noise-free clone realisations (black line; Quartiles). We compare 
this signal with the average signal and its error which is obtained multiply¬ 
ing the standard deviation by 0.67 (red line; Gaussian). In the case of a 
Gaussian distribution, the curves should coincide and the error bars should 
be the same. Bottom panel', same results as top panel but using the noisy 
version of the clone. 


ing to the field of view of each simulation), given the particular set 
of cosmological parameters in the clone. 

The distributions are shown in Figure [7] for angular scales 
0 = 3.4,30.0 arcmin of the equilateral The functional 

form of these distributions depends on the sampling and the pres¬ 
ence of shape noise (noisy). In order to better illustrate the non- 
Gaussian feature of the likelihoods we plot their best-fitting Gaus¬ 
sian inside the panel; the Gaussian fits have the same mean and 
variance as our observed values. In this figure, we find clear devi¬ 
ations from a Gaussian model in the absence of shape noise, the 
observed distributions are skewed towards large values of (Mf p ) 
(bottom panels). The skewness becomes weaker in the presence of 
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Figure 6. Measurements of the EEE (Mf p (0i, 02 , 63 )) in the 120 CFHTLenS pass fields (dark green diamonds). The aperture radii 0 1: are quoted in units 
of arcmin. The error bars have been computed using 184 sight lines of the clone. The average signal in the simulations is shown in red (solid lines or filled 
circles), and the analytical WMAP5 predictions with the SC01 bispectrum are shown in blue (dashed lines or stars). Each panel show the EEE signal as a 
function of 63 62 for a fixed value of 0\ and 9 2 . Panels on the diagonal correspond to equilateral configurations with, increasing size from top to bottom. In 

each row, from left to right we increase the size of 9 2 , keeping 6 1 fixed. The filled orange boxes are the best-fit (fl m , erg) = (0.32, 0.7) of the SC0I model 
to the CFHTLenS data for 9i ^ 5 arcmin. For details see Sect. [ 6 ] 


shape noise but deviations from the Gaussian description are still 
discernible, especially for 9 = 30 arcmin (top right panel). 

For the range of angular scales considered in this paper, the 
non-Gaussian features are most prominent for 9 > 10 arcmin as 
highlighted by Fig. [8] In this figure, we plot for Gaussian the av¬ 
erage and for Quartiles the median of (Mf p (0)) as function of 9 
using the simulations. In addition, we indicate as error bars the stan¬ 
dard deviation times 0.67 (Gaussian) and the ±25 percent quar¬ 
tiles ( Quartiles ). In the case of a Gaussian distribution of (Mf p V 
the data points and error bars of Gaussian and Quartiles should be 


identical, whereas differences indicate deviations from a Gaussian 
distribution for a given 9. Such deviations become most obvious 
for 9 > 10 arcmin, which could indicate that we have to find a 
non-Gaussian model of the likelihood for a cosmological analysis 
of CFHTLenS. To test this hypothesis in the following, we utilise 
the clone distribution of (Mf p ) of the noisy sample to constmct a 
numerically sampled non-Gaussian likelihood and compare its con¬ 
straints on (fl m , ag) to the constraints from a traditional analytical 
Gaussian model. 
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5 COSMOLOGICAL ANALYSIS 

In this section, we describe the details of the cosmological analy¬ 
sis of (Map). For the analysis we use Nd = 20 combinations of 
aperture radii between 5 < 8 i ^ 30 arcmin in (M® p (0i, 62 , 83)), 
shown in Figure [6] We combine all Nd data points into one data 
vector d. The key features of our analysis are a non-Gaussian like¬ 
lihood function that we estimate from the clone and a data com¬ 
pression to suppress the numerical noise in the final result. We com¬ 
pare the results of this analysis to the results with a standard Gaus¬ 
sian likelihood. Moreover, our new approach allows us to factor in 
model uncertainties as well as measurement noise. We include er¬ 
rors due to shape noise, cosmic variance, measurement errors of 
source ellipticities, uncertainties in the Hubble parameter, and a 
multiplicative error in the overall amplitude of the predicted matter 
bispectrum. We embed everything into a Monte-Carlo scheme for 
which no analytical form of the likelihood has to be specified; only 
realisations of statistical errors have to be provided. 

5.1 Model parameters and priors 

For the comparison of our non-Gaussian likelihood to a standard 
Gaussian likelihood, we jointly constrain the two cosmological pa¬ 
rameters of a flat ACDM cosmology with a cosmological constant 
F2 a = 1 — fl m : the matter density and the r.m.s. dispersion og 
of matter fluctuations on a scale of 8/i -1 Mpc linearly evolved to 
the Universe of today. Moreover, we assume for the Hubble param¬ 
eter Ho = h 100 kms -1 a Gaussian prior with mean h = 0.7 and 
relative variance of 5%, compatible with the CMB-only constraints 
of WMAP5 for a flat universe. We compile the model parameters 
into the parameter vector p = (fl m , cr g,h). All other cosmological 
parameters are fixed to their WMAP5 best-fit values to be consis¬ 
tent with the clone; see Sect. 13.II and 13.21 for all relevant model 
parameters. We denote all priors by the probability density function 
(PDF) P„(p). 

In addition, we assume that the overall amplitude of (M® p ) 
can only be predicted up to multiplicative factor 1 + / with a Gaus¬ 
sian prior Pf(f) of r.m.s. 20% and mean (/) = 0. This directly 
accounts for the spread in predictions for the matter bispectrum by 
competing models (see Fig.[3j. 

5.2 Likelihood 

Let d be a vector of observables in an experiment, such as our mea¬ 
surements of (M|p). For the statistical analysis of these data d, we 
can distinguish between two categories of errors. 

Firstly, errors originating from the theory side, or theory er¬ 
rors , which depend on the model parameters p in general. Theory 
errors are present if a model m(jp) cannot perfectly fit d in the to¬ 
tal absence of measurement noise. In our analysis, three kinds of 
theory error are accounted for: cosmic variance, intrinsic shapes of 
the sources, and model uncertainties in the (cosmic average) matter 
bispectrum. Cosmic variance arises because our model predicts the 
cosmic average of (Mf p ) but not the statistics for a particular vol¬ 
ume of the Universe. Intrinsic shapes of sources are not observable 
although they can practically be inferred with reasonable assump¬ 
tions. Hence intrinsic galaxy shapes are nuisance parameters in our 
model that we marginalise over. And, as reflected by Fig. [3] there 
is a set of theoretical models for the matter bispectrum currently 
available that vary among themselves by about 20% in amplitude 
for the range of angular scales that we study here. We parametrise 
this particular theory uncertainty by the nuisance parameter /. 


Secondly, even if a model fits the data we expect residuals 
of the fit due to errors in the measurement process: the measure¬ 
ment noise. Measurement noise is by definition only related to the 
observables and hence unrelated to p. Our main source of measure¬ 
ment errors, and the only one we account for here, are uncertain¬ 
ties in the galaxy ellipticities et. Other conceivable sources such as 
source positions or their redshifts are neglected here. 

In our model of d, we include theory and measurement errors 
by the Ansatz 

d = (1 + f)m(p) + A d(p) =: m f (p) + A d(p) , (39) 

where / is the multiplicative error in the model amplitude m(p), 
while A d(p) comprises all remaining theory and measurement er¬ 
rors combined. In more complex applications, theory errors and 
measurement errors could be separated and simulated indepen¬ 
dently from each other. On average we have (A d(p)) = 0 and 
(/) = 0. Our model of the likelihood function is then as follows. 
We denote by P/sd{Ad\p) the probability density of an error Ad 
given the parameters p. To include the error of /, we then write the 
likelihood C(d\p) of d given p as 

C(d\p)=J df Pf(f)PAd (d — m/(p) . (40) 

Up to a normalisation constant E(d) the posterior PDF of (H m , og ) 
given d is therefore 

Pp(H m , cr 8 |d) = E^ 1 (d) J dhP p (p)£(d\p) . (41) 

The exact value of the so-called evidence E(d) is irrelevant for this 
analysis and hence set to unity. The Hubble parameter h is weakly 
constrained by the data. Therefore, we marginalise the posterior 
over h. 

A reasonable and common approximation of the likelihood 
C{d\p) would be a multivariate Gaussian PDF as, for instance, 
applied in FI 4. As shown in Sect. 14.31 by our simulation of the 
CFHTLenS data set, however, the (M® p ) measurement exhibits 
deviations from Gaussian statistics. This motivates us to apply a 
non-Gaussian model of the likelihood. Towards this goal, we out¬ 
line in the following an algorithm that estimates £(d\p) based on a 
discrete set of Monte-Carlo realisations of A d(p) for the posterior 
Pp(H m , (7g|cl). 

5.3 Monte-Carlo sampling of likelihood 

An opportunity for an approximation of C(d\p) by a Monte-Carlo 
process can be seen by re-writing the PDF P&d as 

P/sd(Ad\p) = J dAx P/sd(Ax\p)S D (^Ad — Ax'j (42) 

I JVAd / \ 

~ \ Ad ~ Axk ) (43 ) 

k=1 

with 5y>(x) being the Dirac delta function. The sum in last line 
approximates the integral in the first line b y a nume ri cal M onte- 
Carlo integration with AAd points Axk JPress et aill 19921 ). For 
this approximation, we produce, for a given p, N&d random re¬ 
alisations of Ax from P/\d{Ax\p) by means of the clone. We 
denote this process simply by Axt -tv PAd(Ax\p) in the follow¬ 
ing. In the limit N&d —i► 00 , the number density of the points Ax*. 
at Ad converges to Prsd{Ad\p) up to a normalisation constant. 
Even for finite N/s.d, the number density of points still provides an 
useful estimator for the likelihood. This is the basis of our Monte- 
Carlo scheme. For this scheme, we address the infinite noise in the 
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previous estimator by smoothing the number density of sampling 
points. 


5.4 ICA-based interpolation 


Let Aa;fe -r\ -PAd(Aa;|p) be a sample of N Ad simulated dis¬ 
crete points with average (A Xk) =0. In addition, let the ma¬ 
trix N = AatfcAxfc/AAd be an estimator of the covariance 
of Axk- Our aim is to use the smoothed local number density of 
the set of points {Axk} at Ad as approximation of PAd(Ad\p). 
To this end, we exploit the independent component analysis (ICA) 
to factorise the sparsely sampled higher-dimensional likelihood 
P Ad (Ad\p) int o a pro duct of one dimensional (ID) histograms as 


Hartlap et al.l l2009t) . We refer the reader to this paper for a gen¬ 


eral discussion of the ICA and only briefly summarise the practical 
steps here. 

The ICA defines a linear transformation M and new co¬ 
ordinates c := MN _1 / 2 Ad that allow us to express the Ad- 
dimensional PDF 

N d 

P Ad (Ad\p) « detMN ~ 1/2 x P c (i) ([MN~ 1 / 2 Ad]i|p) (44) 


as a product of the ID probability densities P^\ci\p). The con¬ 
stant pre-factor on the right-hand side is the Jacobian of the map¬ 
ping MN^ 1,/2 . Therefore, by means of the ICA all components a 
of c are made statistically independent. Note that for a Gaussian 
P Ad (Ad\p) the ICA is essentially a principal component analysis. 
We estimate the transformation matrix M from {A xk I by applica¬ 
tion of the fast ICA algorithm IE varine nl ll999h . The matrix M 
is a Nd x Aa square matrix in our case, i.e., the number of com¬ 
ponents of c and Ad are the same. 

For the sample {Axk}, the ICA transformation then gives us 
a new set of N Ad sampling points Ck := MN -1 ' 2 Axk- We denote 
the i-components of Ck by Cki = [ Ck]i . For a fixed i, we adaptively 
smooth the ID frequency distribution of these points {cm} to pro¬ 
duce the histogram Q d l A {x\ cm) and use its value at x = d as an 
estimator of P ( a |p). Specifically, to obtain Q[^ a (x ; cm) we ap¬ 
ply a 10th neighbour adaptive smoothing: for a given value x we 
find the distance dio,i(a;) := \x — cio,i| to the 10th nearest neigh¬ 
bour sampling point cio,i € {cm} for fixed i and then compute 

Q ^ {X ’ Ck ' )= d^) ■ (45) 

In doing so, we bin for every index i the ID distribution of values 
{cki} with an adaptive bin width di o.i(x) that depends_on the lo¬ 
cal density of the points {cm} at x dLoftsgaarden & Ouesenberrvl 
Il965l) . We find that this technique is more robust compar ed to 
the Gaussian kernel method with fixed kernel size in lHartlap et all 
d2009h . if the PDF-sampling has a few extreme outliers in the tail 
of the distribution. 

Finally, our approximation of the logarithm of P Ad (Ad\p) is 
given by 


In P A d(Ad\p) rs lnQ ica (Ad; Aaifc) (46) 

n a 

:= 7? + ^lnQ«([MN- 1 / 2 Ad] i ;c, i ) 

i= 1 

where P 1 is the normalisation of Qi Ca (Ad; Ax *,). The normalisa¬ 
tion 7 ) can be ignored in our case because for this study it is identical 
for every position in the parameter space. 


0-1 =5.9 arcmin 10.2 17.5 30.0 




Figure 9. Illustration of the Karhunen-Loeve data compression. All mod¬ 
els assume h = 0.7. Top panel: The red sticks show simulated data 
\ M 2 p ) with Nd = 20 elements sorted by increasing 0\ ^ (I 2 7 If 3 ; 
9i = 5.9,10.2,17.5, 30.0 arcmin. The KL modes are averages of all data 
points with relative weights as indicated by the lines for the mode indices 
n = 1, 2, 3 (right hand j/-axis). Bottom panel: Amplitudes of KL modes 
as function of the mode index for four different noise-free model vectors 
(lines); the open square data point at n = 1 is off the chart at an amplitude 
of approximately -30. The open diamonds show a one noisy clone data 
vector in comparison. 


5.5 Data compression 

The ICA interpolation technique is prone to a bias that under¬ 
estimate^ the size of the credible intervals of the model parameters 
dl I art I an et al.ll2009l) . This bias is large if Nd/N Ad ~ 1 and small 
for Nd/N Ad <C 1; here N Ad — 200 is our number of sampling 
points of the likelihood and Ad = 20 is the size of the uncom¬ 
pressed data vector. A similar problem affects traditional Gaussian 
likelihood analyse s when the inverse covariance has been estimated 
from simulations dHartlan et alJl2007l : lTavlor et alj|2013l) . In order 
to reduce this bias, we compress the size N d of the data vector 
d, and we apply the same compression to the realisations Aa;*, in 
the ICA interpolation. This data compression also has the benefit 
to produce a smoother posterior of p as the sampling noise in the 
interpolated likelihood is reduced. 

Our data compression is based on a Karhunen-Loeve 
(KL) transform iTegmark et alJll997l : iKilbinger & MunshU 120061 : 
lAsgari & Schneideill2014l) . It identifies Ad signal-to-noise (SN) 


© 0000 RAS, MNRAS 000, 000-000 






































































CFHTLenS: non-Gciussian analysis of 3pt shear correlations 13 


modes of d and rejects modes with a low SN. This procedure al¬ 
lows us to find a linear projection d' = C d that reduces the number 
of elements of d to N# < Nd. The projection matrix C is chosen 
such that most of the cosmological information on p is still avail¬ 
able in d'. Importantly, the KL transform - or any projection ma¬ 
trix C for that matter - at most decreases the constraints on p but 
does not bias the result provided the data can be fit by the model. 
The KL transform affects the cosmological analysis and the ICA 
decomposition only in that far that we now use as data vector d' 
instead of d, the compressed model vector m!(p) = C m(p), and 
the Monte-Carlo realisations Ax' k = CAx k . 

For a data compression matrix C, we would like to identify lin¬ 
ear combinations of the components of d that vary strongly when 
changing the parameters p (signal variance). On the other hand, 
in the presence of measurement noise and theory errors, we would 
also like to identify combinations that are most significant com¬ 
pared to noise. The KL transform finds a compromise of both. This 
means the KL transform identifies modes ei in d-space that have 
large ratios of signal variance and noise variance. The compression 
matrix C = (ei,..., ejv d , ) T comprises the Nd' KL modes with 
the highest SN ratio. To construct C, we proceed in two steps. 

First, we compute the signal covariance S of model vectors 
m(p) over a volume V in parameter space that is defined by our 
prior information P P {p) on the cosmological parameters, 


f dQ^dasdh P p (p) , , w , , 

S= I - y - (m(p) - (m))(m(p) 

where V = f dfl m dasdhP p (p), and 

, , f d£l m da$dh P p (p) 

(m) = / ---— -m(p) 


( m )) T , 

(47) 

(48) 


is the mean of m{p). For the computation of S by integration, we 
adopt a flat prior P p (p) for € [0.1, 1.0], as € [0.4, 1.0], and 
h £ [0.6, 0.8]. The previous integrals could in principle be esti¬ 
mated by a Monte-Carlo process for which we (i) randomly and 
repeatedly draw parameters p from the prior P p (p ), (ii) compute 
m(p), and (iii) determine the mean (m) and covariance S of all 
realisations of m(p). The Eigenvectors of S with the largest Eigen¬ 
value determine modes of d that are most sensitive with regard to 
P- 

Second, we compute the noise covariance, which is simply 
the covariance matrix N of Ax for a chosen fiducial cosmology 
(Sect. lL4l. The data compression is thus optimal for the fiducial 
model. The KL modes are then the Eigenvectors e, of the gener¬ 
alised Eigenproblem 


Sei = AiNei, 


(49) 


where A i are the Eigenvalues of e t ; the SN ratio of the mode e t is 

x/A i. 

In the top panel of Fig. [9] we show simulated data (M^ p \ for 
our Nd = 20 combinations of aperture radii (red bars). The KL 
modes are linear combinations of the aperture radii with relative 
weights defined by the components of et. These weights are dis¬ 
played as lines for the first three modes (right hand y-axis). Im¬ 
portantly, the weights can have positive and negative signs. The 
lines show that our KL modes are mostly sensitive to large aperture 
scale radii with 0i ^ 10.2 arcmin. The bottom panel shows as lines 
four compressed model vectors Cm(p) that do not contain noise 
in comparison to one noisy clone data vector. The error bars re¬ 
flect the la uncertainty for noise as in our CFHTLenS data. Except 
for models with large values of f 2 m and as (grey dashed-dotted 
line) all models are essentially zero for KL mode indices larger 


than three. This means their cosmological information is mostly 
inside the first couple of modes. Note that errors between differ¬ 
ent KL modes are uncorrelated and of equal variance by construc¬ 
tion, i.e., (Ad'Ad' ) = <5]^. Nevertheless, in the case of a non- 
Gaussian statistics, there may be higher-order correlations such as 
(Ad'iAd'jAd'ff ^ 0 present in the data. 


5.6 Bias correction 

After a data compression with N& ^10 and the following correc¬ 
tion of the likelihood, we expect the bias in the size of the credible 
interva ls of our samp l ed po sterior to be negligible for the following 
reason. lHartlap et al] ( l2009l) report that this bias is similar to the one 
known for a Gaussian likelihood for which the data covariance is 
estimated from realisations of the data. In the Gaussian case, to de¬ 
bias the likelihood the inverse of the estimated covariance has to be 
rescaled by a factor/ corr = (AAd — TVa/—2)/(7VAd — 1 ) or, equiv¬ 
alently , L LAorr \y e a pp]y the latter correction to our sampled 
likelihood Qi ca ,(Ad'; Ax k ) in Eq. 146L However, this correction 
is small since f COTT ^ 0.94 for Nd' ^ 10. 


5.7 Sampling algorithm of posterior 

We compute the posterior P p (fl m , crs|rf) on a 100 x 100 grid that 
covers the (fi m , <7s)-plane. To compute the posterior value for a 
grid pixel we proceed in four steps: 

(i) Generation of a set of N A d realisations d k to sample the error 
P Ad (Ad\p) for a given p: We obtain the d k by measuring (Mf p ) 
on the noisy-sample of the clone, hence both the intrinsic galaxy 
shapes and the cosmic variance are simulated. In the simulation, 
the ID variance of intrinsic source ellipticities is set to <r e =0.28 
to account for both the true intrinsic variance and the measurement 
error of ellipticities; both are assumed to be Gaussian. From the 
clone patches, we compute a simulated (M^p) for an effective 
CFHTLenS survey area of roughly 100 square degree by combining 
the clone measurements of seven randomly selected sight lines. 
Each clone line-of-sight, out of 184 in total, covers 12.84 deg 2 . 
Repeating this procedure N A d = 200 times gives us the sample 
Ax k = d k — d where d expresses the average of all realisations 

d k . 

(ii) Determination of the data compression matrix C for a given 
covariance N of {Acc*,}, models m(p), and the model prior 
Pp{p)'. see Sect. 15.51 for details. From this we compute the com¬ 
pressed vectors d! = C d and Ax' k = CAx k . We compute the 
compression matrix C once for a fiducial cosmology and use it for 
all p. 

(iii) ICA factorisation of the compressed P Ar i' (Ad'\p) into a 
product of interpolated, ID histograms Qffx; cm) by means of 
an ICA of the sample c k i = [Ax' k \i\ see Sect. l5.4l for details. 

(iv) Marginalisation over the multiplicative factor / and the 
Hubble parameter h to obtain the posterior value P p (Q. m , as\d'): 
We perform the marginalisation by another Monte-Carlo integra¬ 
tion that draws IVfh = 500 values fi Pf(f) and the same 
amount of values hi N( 0.7, <Th) from a Gaussian prior for the 
Hubble parameter. With these values we compute 

1 N lh 

P p (Q m ,a s \d') » [Qi C aK-mUPi);Aa4)] /o ° rr , 

i= 1 

(50) 

where m'f i (p i ) = m'(jPi)(1 + 0.2/*) is the rescaled 
model m , (p i ) for p i = (fi m , as, hi). The approximation 
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Qi C a(A<f'; Ax' k ) depends on the factors Q[^{x] Cki), see Eq. 

Ideally, step (i) is repeated for any new value of p based on 
a simulation for a fiducial cosmology with parameters p. To keep 
the simulation effort viable, however, we assume that the scatter 
of dk is as in the clone for any p in the entire parameter space 
explored. Therefore step (i) only has to be performed once. For a 
Gaussian likelihood, this assumption would be equivalent to using 
the same likelihood covariance for all p, which is in deed a com mon 
ass umption as e.g. in F1 4 (for a discussion see lEifler et all 120091 
and lKilbinger et ai]|2013ll . In a future application of our technique, 
subject to less prohibitive computational constraints, a computation 
of the sampled likelihood for different p is conceivable. 

Finally, it is likely that we moderately underestimate the scat¬ 
ter of Axk in step (i) due to cosmic variance. The 184 realisa¬ 
tions of the clone allow only for 26 independent realisations of a 
100 deg 2 survey; the cosmic variance between the Axk is partially 
correlated. 


5.8 Tests of the posterior construction 

In this section, we test the robustness and accuracy of our code 
implementation of Sect. l5.7l before its application to the CFHTLenS 
data. 

We start by looking at different degrees of data compression 
in the top panel of Fig. EH For this purpose, we use a noise-free 
model vector m(p) with the parameters p = = 

(0.3, 0.8, 0.7) as input to the cosmological analysis. We plot three 
credible regions of the resulting posterior P p (G m , ag\d). The max¬ 
imum of all posteriors coincides with the parameters of the input 
vector (not shown). The figure supports our earlier expectation that 
only few KL modes are required for the analysis: there is no clear 
improvement for N,y ^ 3. The differences between the panels are 
likely due to numerical noise in the sampled likelihood for the fol¬ 
lowing reason. In Fig. m we plot the posterior of the combined 
quantity Es = <J8(fI m /0.27) 0 ' 64 for the previous model data vec¬ 
tor but devising now a Gaussian likelihood. This likelihood model 
is not subject to sampling noise as its analytical form is completely 
determined. As can be seen, the posterior P P CEs\d') is virtually 
unchanged for N& ^ 3. Therefore, the first two modes contain all 
cosmological information of the model data. The same is presum¬ 
ably true for the non-Gaussian case so that all variation in the top 
panel can be attributed to numerical noise. 

As can be seen in Fig. [9] (unlikely) models with values high in 
both and as exhibit some signal until the mode n = 7. There¬ 
fore, we set N^i = 5 in our final CFHTLenS analysis as compro¬ 
mise between sampling noise and cosmological information. More¬ 
over, the middle panel shows that the inclusion of a 20% systematic 
error in the model amplitude has a small impact on the posterior 
information (solid black line). It slightly stretches the uncertainty 
toward larger values of Es. 

For the bottom panel of Fig. EH we test the accuracy of our 
analytical matter bispectrum and that of the likelihood analysis by 
means of the clone. We run our analysis code for 26 independent 
simulated CFHTLenS-like data vectors created using noise-free re¬ 
alisations; we combine the 184 patches into 26 groups of seven data 
vectors. The figure shows three credible regions of the combined 
posterior of all 26 analyses. We attain the combined posterior by 
adding the 26 individual logarithmic posteriors on a grid. Ideally, 
the resulting posterior should be consistent with the cosmological 
parameters in the clone that we indicate by the intersection of the 
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Figure 10. Results of our code tests running on mock data. Top panel: A se¬ 
ries of cosmological analyses using different degrees of data compression. 
All four runs use same simulated data vector with = 0.3, erg = 0.8, 
and h = 0.7 but a different number of KL modes as indicated. The coloured 
regions highlight the 68%, 95%, and 99% credible regions for a flat cos¬ 
mology. Insets with more KL modes N^/ are subject to more numerical 
noise. Middle panel: The posterior of Es = <78(H m /0.27) 0 ' 64 for the test 
data vector for different degrees of compression and a Gaussian likelihood. 
The dotted black solid line uses N&r = 5 and no error in the model am¬ 
plitude (/ = 0). Bottom panel: Results of a verification run of the anal¬ 
ysis code. The credible regions reflect the combined constraints from 26 
simulated noise-free CFHTLenS measurements in the clone simulation. 
The posterior is offset with respect to the fiducial model (Q m = 0.279, 
erg = 0.817; dashed black lines). 
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two dashed lines. The contours, however, indicate with about 3cr a 
model bias due to which we do not perfectly recover the clone 
cosmology: the centre of the contours, indicated by the black point, 
is offset in <78 by a few percent and in by roughly 10 percent. 
However, this offset is still small compared to the CFHTLenS noise 
levels, which can be seen by comparing offset in the bottom panel 
to the size of the credible regions in the top panel. To assure that 
the bias is unrelated to either the data compression or the ICA inter¬ 
polation of the likelihood, we also determine the maximum of the 
posterior for the average, uncompressed clone data vector for a 
simple Gaussian likelihood. We find biased values consistent with 
the non-Gaussian results of compressed data. Our interpretation is 
that the bias originates from a small mismatch between the analytic 
model and the clone average that is still visible at 5 arcmin (EEE 
mode in Fig. [2}. 

In conclusion, there is a systematic bias in our analytic model 
for the third-order aperture statistics as revealed by our comparison 
to N-body simulated data. Since this bias is small relative to the 
levels of shape noise and cosmic variance in CFHTLenS we can 
ignore it for the scope of this paper. For the analysis of next gener¬ 
ation surveys, on the other hand, more accurate models of (M® p ) 
will be needed. 


6 RESULTS 

We present our constraints of ag and for the CFHTLenS 

data for a flat ACDM model in Fig. QT] Therein we include 
only measurements of the EEE mode of 1 , 62 , 83 )) for 

9i = 5.9,10.2,17.5, 30.0 arcmin. In order to avoid confirmation 
bias in the generation of this result, we blindly analysed four data 
sets simultaneously of which three were noisy mock data and one 
was the CFHTLenS data vector with the results shown here. The 
CFHTLenS data vector was revealed only after the cosmological 
analysis. Furthermore, we performed the analysis after the system- 
atics checks of the data and after our final decision on the usable 
range of angular scales. In particular, we did not use the CFHTLenS 
data during the development and testing phase of the analysis code 
in any way (Sect. [5}. 

In the top panel of Fig. M we show the 68%, 95%, and 
99% credible regions for the joint constraints of the matter den¬ 
sity parameter and the amplitude os of the matter density 
fluctuations. We performed the analysis twice: once with the non- 
Gaussian model of the likelihood for the left inset (non-Gaussian), 
and once with a Gaussian likelihood for the right inset (Gaus¬ 
sian). The Gaussian likelihood is based on the noise covariance 
matrix N in Sect. ns Both posteriors use the same degree of 
data compression (five KL modes). For reference we have in¬ 
dicated the slightly differing best-fitting value s of WMAP9 and 
Planck as solid and dashed lines, respectively fainshaw et alj20l3 : 
IPlanck Collaboration et al.ll2014l) . In summary, there is little dif¬ 
ference in the posteriors, mostly visible between the 68% credible 
regions. Both constraints are highly degenerate and basically only 
exclude simultaneously large values of S2 m and ag¬ 
in order to break the degeneracy of the parameters for the 
bottom panel of Fig. m we assume an additional prior on ei¬ 
ther (left) or os (right). The top-hat priors are centred around 
fl m = 0.275 and <7s = 0.8, close to the best-fitting values of the 
WMAP9 results. The widths of the top-hats are Afi m = 0.05 
and Acrs = 0.1. The resulting marginalised posteriors of the 
non-Gaussian data model are shown as grey bars. With these 
strong priors we infer = 0.27(tg'o5 and <7s = (68% 




Figure 11. Top panel: Constraints in the <78 — fl m plane from the 
CFHTLenS (Mf p ) for a flat ACDM cosmology with 5% uncertainty on 
h = 0.7 and a 20% uncertainty in the model amplitude (credibility red: 
68%, orange: 95%, yellow: 99%). The dashed lines indicate the Planck 
best-fit parameters, the solid lines WMAP9 parameters. The left inset de¬ 
picts the posterior for a non-Gaussian likelihood, whereas the right inset is 
based on a Gaussian model. The dotted line on the left (right) corresponds 
to the best-fitting value of Yg = ag (%,„/().27)°' f ' 4 = 0.8(0.73). Bottom 
panel: Marginalised cosmological results for ag (left) and Q In (right) alone 
with different top-hat priors on the parameters (as indicated). The thick blue 
dashed line indicates constraints from the Gaussian likelihood in compari¬ 
son to the non-Gaussian likelihood (grey area). The thick dotted black line 
are non-Gaussian constraints for a fixed h = 0.7. whereas the thick dotted- 
dashed green lines / = 0 do not include the 20% multiplicative error in the 
model amplitude. Except for the red dashed lines, the posteriors are nor¬ 
malised to their value at the maximum. 


credibility). Furthermore, without the imposed top-hat priors the 
marginalised posterior of either parameter yields only weak con¬ 
straints (thin dashed red lines), which are essentially just the up¬ 
per limits < 0.45(0.67) and<7 8 < 0.75(0.93) for a 68% (95%) 
credibility. These upper limits depend on the adopted broad priors 
of 0.15 ^ ^ 1 and 0.4 ^ ag 1 due to the strong degener¬ 

acy of both parameters. Therefore, it is more relevant to combine 
(H m , <7 8 ) into E 8 = <78(f2 m /0.27) 0 ' 64 as done in F14. We plot the 
posterior of E 8 for the non-Gaussian data model in the bottom 
panel of Fig. ( ITTt and obtain E 8 = 0.79lto!ii (grey area). 

The constraining power of the non-Gaussian model for low 
values of cr 8 < 0.7, < 0.2, and S 8 < 0.7 is slightly better 

compared to the Gaussian data model, as can be seen in the bot¬ 
tom panels of Fig. [TO and Fig. |T2] (compare the grey histogram 
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Figure 12. Top panel: The first five KL modes of the compressed 
CFHTLenS data vector (red triangles) in comparison to both the best-fitting 
model (filled blue diamonds) and the scatter of amplitudes of our mod¬ 
els for a parameter range (grey whisker boxes). For the latter, the filled 
boxes denote the amplitude spread of 68% of all models around the median 
for Eg = cT8(O m /0.27) 0 '® 4 ^ 1.0 and h = 0.7, whereas the whiskers 
bracket the entire amplitude range. There is a strong tension between mod¬ 
els and data for n = 3. The error bars do not include the 20% multiplicative 
error in the model amplitude. The \ 2 per degree of freedom of CFHTLenS 
is 4.2 or 2.9 when the multiplicative error is included. The open diamonds 
denoted clone correspond to one random realisation of a noisy clone mea¬ 
surement. Bottom panel: Posterior of Es for the non-Gaussian data model 
( non-Gauss ) for various strong priors (lines / = 0 and h = 0.7) compared 
to the posterior in the Gaussian model (line Gauss). The red dotted-dashed 
line shows the posterior based on the KL modes n = 1, 2, 4, 5 only. 

to the blue dashed line). In addition to that, the black dotted lines 
shows our non-Gaussian constraints for a fixed h = 0.7 without 
the marginalisation of the Hubble parameter that all other results 
are subject to. The difference to the posterior in grey is small. Like¬ 
wise, the impact of the 20% error in the predicted model amplitude 
is also small, as indicated by the dotted-dashed green lines (/ = 0). 

Figure [T2] displays the goodness of the SC01 model with re¬ 
spect to the CFHTLenS data. We plot the first five uncorrelated KL 
modes that are used in our analysis (red triangles). For comparison 
we also plot the KL modes of one random noisy clone data vec¬ 
tor into the plot (open diamonds). The blue filled diamonds corre¬ 
spond to the best-fitting model of CFHTLenS with = 0.32 and 
os = 0.7. We observe a clear discrepancy between the best fit and 
the CFHTLenS data at n = 3 and some weak discrepancy at n = 2. 
The \ 2 °f the residuals of the best fit is 4.2 per degree of freedom 


or 4.2/(1 + 20%) 2 = 2.9 if we account for a 20% systematic error 
in the predicted amplitude; the systematic error can be included by 
increasing the errors of modes in the plot by 20%. In order to illus¬ 
trate that reasonable models are unable to fit the CFHTLenS data at 
n = 3, we have added the grey boxes to this figure. They depict the 
scatter of 68% of the model amplitudes for parameters Es ^ 1.0 
and (fl m , as) £ [0.1,1] x [0.4, 1]; current data constraints confine 
Eg to about 0.76 ± 0.06 with 68% confidence (F14). Additionally, 
the whiskers bracket the total amplitude range covered by the mod¬ 
els. Clearly, our n = 3 mode is well outside the whisker region, 
underlining no model reproduces our measurement here. Note that 
mixed higher-order moments between the errors of the KL modes 
may exist even though their second-order correlations vanish. 

Finally, the red dotted-dashed line in the bottom Fig. [12] cor¬ 
responds to the posterior of Es leaving out the KL mode n = 3. 
Statistically the posterior is consistent with the constraint for non- 
Gauss although the mode of the distribution shifts to slightly 
smaller values of Es and tightens a bit, i.e., Es = 0.77 ± 0.08. 


7 DISCUSSION 

We conclude from our results that the model of a Gaussian likeli¬ 
hood is sufficient for the cosmological analysis of the CFHTLenS 
{Ml p) data, or data of comparable surveys, despite the evidence 
for a non-Gaussian distribution of errors. To arrive at this conclu¬ 
sion we measured the third moment (Mf p ) of the aperture mass 
for aperture scales between 5 and 30 arcmin. For this we excluded 
scales below 5 arcmin owing to indications of systematics in the 
aperture statistics. Our choice of 5 arcmin as smallest usable scale 
in the analysis had been made prior to the cosmological analysis 
of the data to avoid confirmation bias. Our maximum scale is de¬ 
termined by the geometry of the survey. For this angular range, 
the clone simulation of our data in Fig. [7] and in Fig. [^indi¬ 
cates a moderate non-Gaussian distribution of the measurement 
errors between 10 and 30 arcmin. A Gaussian likelihood hence 
may produce biased parameter estimates in a cosmological anal¬ 
ysis. To test the validity of a Gaussian likelihood, we performed 
a comparative cosmological analysis with a non-Gaussian likeli¬ 
hood that is constructed from the distribution of (Mf p ) estimates 
in our simulation. For this task, we devised a novel technique that 
involves a data compression to reduce the aperture statistics to a 
few essential uncorrelated modes that are shown in Fig. [T2] (higher- 
order correlations may be non-vanishing). Our statistical analysis 
factors in both measurement noise and uncertainties from the the¬ 
ory side without the need to specify the analytic form of the error 
PDF. This practical algorithm is an important contribution of this 
paper for future analyses because its applicability is not restricted 
to third-order shear statistics. The outcome of our statistical anal¬ 
ysis for the two cosmological parameters (fi m , crs) is displayed in 
Fig. II Hand juxtaposed with a second analysis that uses a Gaussian 
likelihood. The posteriors of the parameters exhibit only little dif¬ 
ferences between the Gaussian and non-Gaussian likelihood - at 
least for a flat AC DM model which i s stro ngl y fav oured b y recen t 
CMB constraints IPlanck Collaboration et alJl20l4 rHinshaw et al. 
120131) . The main difference between both models of the likelihood 
is illustrated by the bottom panels of Fig.QT|and Fig.[l2] the non- 
Gaussian model excludes more strongly low values of as, the lower 
68% bound shifts up from 0.63 in the Gaussian case to 0.66 in the 
non-Gaussian; the lower limit of Es shifts from 0.60 to 0.68; the 
improvement of the lower limit of is below 0.01. A qualitative 
similar behaviour for and as was found with the non-Gaussian 
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model of lHartlap et al.! j200gh for the second-order statistics of cos¬ 
mic shear. The observed changes are small compared to the width 
of the posteriors though. Consequently our results support the re¬ 
cent decision of F14 to apply a Gaussian model of the likelihood to 
angular scales smaller than 15 arcmin as reasonable approximation, 
although their cosmological constraints can probably be tightened 
by a non-Gaussian data model. 

For future surveys with increased survey area and higher 
source number densities, we expect that non-Gaussian features will 
become more prominent in the data; a Gaussian model may hence 
then no longer be adequate. As seen in the top panels of Fig. [7] on 
small angular scales around 5 arcmin the signal is dominated by 
shape noise of the sources, which tends to Gaussianise the distri¬ 
bution of (M^f) for statistically independent intrinsic shapes. On 
larger scales of about 30 arcmin, on the other hand, cosmic variance 
is dominating, which on our sub-degree scales is still non-Gaussian 
due to the non-linear clustering of matter (Fig. [8J. For a fixed an¬ 
gular scale and for an increasingly higher number of source triplets 
in the estimator, the amplitude of cosmic variance grows greater 
relative to the shot noise variance. Flence future lensing surveys 
will exhibit stronger non-Gaussianities on the angular scales that 
are considered here. However, for a quantitative assessment of the 
impact of these non-Gaussianities on parameter constraints with 
~ 10 3 square degree surveys we need more independent simulated 
data vectors than the clone can provide. The clone allows only 
for about two independent realisations of a (Mf p ) measurement in 
this case. This is far too few to sample the data likelihood. Still, 
considering that we already see an inaccuracy with a Gaussian data 
model in CFHTLenS, we anticipate that a Gaussian model will po¬ 
tentially bias the cosmological analysis in future surveys. 

Our posterior for (O m ,ag) from (Mf p ) only is highly de¬ 
generate, and it is 68% consistent with recent constraints from the 
CMB and the Gaussian CFHTLenS analysis of F14. The degener¬ 
acy is shown in Fig. eh The joint constraints on two parameters 
naturally have to be degenerate along a one-dimensional line, as 
the cosmological information in the data essentially consists only 
of one significant KL mode, see Fig. |T2] From the joint, degen¬ 
erate constraints we infer Es = crg(f2 m /0.27)°' 64 = 0.797g' < j l j, in 
accordance with F14. Note that F14 used only cosmological scales 
between 5.5 and 15 arcmin as well as a more accurate model for 
the matter bispectrum which explains minor difference between 
our posterior contours. For a comparison to WMAP9 results, we 
break the parameter degeneracy by imposing a narrow prior on 
either parameter or ag that is consistent with WMAP9, i.e., 

£ [0.25,0.30] or as £ [0.75,0.85], For a prior on one pa¬ 
rameter, the constraint of the other parameter is then also consis¬ 
tent with the WMAP9 best fit (erg = 0.82 and f= 0.28): we 
fin d O m = 0.27jln' ni] and ag = 0.77to ii- The recent results from 
IPlanck Collaboration et al.1 ( l2014lt . on the other hand, favour some¬ 
what larger best-fit values a g = 0.83 and = 0.31 due to a 
smaller h = 0.67 ± 0.014. These values are also consistent with 
our CFHTLenS findings although there appears to be a 1 <j ten¬ 
sion between the values of H m . This tension is, however, compa¬ 
rable to the inaccuracy in our statistical methodology or model as 
seen in our verification run (bottom panel of Fig.llOt; and ag 
are somewhat underestimated in the analysis. In conclusion, our 
(Mf p ) cannot distinguish a WMAP9 from a Planck cosmology; 
both are equally consistent with our CFHTLenS results. 

In spite of the broad consistency of (fi m , as) with the stan¬ 
dard cosmology model, the CFHTLenS (Mf p ) data have features 
that cannot be explained by any of our SC01 models or more accu¬ 
rate models of the matter bispectrum in dark-matter-only universes 


(flat ACDM). A dark-matter-only model consequently is not good 
enough to describe our observed third-order correlations in the cos¬ 
mic shear field, or there are remaining systematics in CFHTLenS 
on scales greater than ~ 10 arcmin that are relevant for the third- 
order statistics. The only moderately good model fit to the data be¬ 
comes easily obvious when we inspect the compressed CFHTLenS 
data in Fig. EH error bars are uncorrelated in this representation, 
and they have equal variance; the best-fit has f 2 = 4.2 per degree 
of freedom without 20% error in the model amplitude and f 2 = 2.9 
including the amplitude error; the number of degrees of freedom 
is 3. All models with a generous cut of Eg ^ 1.0 are essentially 
zero for KL modes n ^ 3. The CFHTLenS data at n = 3, on the 
other hand, is 3 a away from any of these models. Only models 
that are simultaneously large in ~ 0.7 and ag ~ 0.8 get closer 
to the n = 3 data point but quickly move away from n = 1 and 
n = 5 at the same time and cannot explain our observation either; 
see e.g. the dashed-dotted line in the bottom panel of Fig. [9] The 
discrepancy is also visible in Fig. [6] where we plot the best-fitting 
model (filled squares) in comparison to the CFHTLenS data (open 
diamonds): the CFHTLenS data points are systematically above 
the model data points for 62 = 5.91 arcmin (third column) in or¬ 
der to be consistent for larger values of 62 (columns 4-5). Cer¬ 
tainly, this disagreement with theory, reflected by the n = 3 KL 
mode, could be a shortfall of SCOT However, no analytic model 
can be better than the bispectrum power in cosmological simula¬ 
tions that are used to either test or calibrate bispectrum models in 
the non-linear regime dValageas et al.ll2012l : iGil-Marin et aljEoTH : 


IScoccimarro & Couchmanl 200lh . Therefore, improved modelling 
would at best reproduce the clone data points which are plotted 
in Fig.[9] open diamonds in the lower panel, and Fig.EH filled trian¬ 
gles in the upper panel. Clearly, these data points do not exhibit the 
n = 3 feature seen in CFHTLenS; we find the same for the clone 
data points in the other 183 line-of-sights. Since both the clone 
and the CFHTLenS two-p oint statistics are consistent with the stan- 
dard cosmological model (Fu et a l .120141 ; Kitching et al] |2014ll . we 
therefore conclude that also a bispectrum model more advanced 
than SC01 or a dark-matter-only standard ACDM in general cannot 
explain our n = 3 mode of CFHTLenS. The n = 3 mode is most 
sensitive at (6»i, d 2 , 63 ) = (5'.9,17'.5, 30') and (10'.2,17'.5,17'.5) 
(top panel of Fig. We hence broadly locate the discrepancy be¬ 
tween 6 ~ 10 — 20 arcmin, or correspondingly t « 248 — 496. 

Despite the evidence of some shear systematics in the data, 
we are unable to conclusively identify either intrinsic alignments 
(IA) or residual PSF systematics as origin of the model discrep¬ 
ancy. However, there is evidence for IA playing only a minor role 
in this context. After the application of the systematics test in Ap- 
pendix[A] we additionally rejected nine CFHTLenS fields from the 
H12 sample. This mildly affects the third-order aperture statistics 
between 10-20 arcmin (left panel of Fig. 0: the pass fields have a 
higher signal for the equilateral We hence suspect that PSF 

systematics are at least in part responsible for a signal deficiency 
around 15 arcmin. Furthermore, after the removal of the early-type 
galaxies from our shear catalogue, we found little difference in 
the EEE signal but rather an increase of the BBB signal on small 
scales (Fig. [5}. This behaviour is the opposite of what is expected 
for IA contaminate d data: most of the IA signal is associated with 
early type galaxies jMandelbaum et al.l2006l ; Ijoachimi et aljEoTH : 
iMandelbaum et ahluOl 1 1. Therefore it is unlikely that the signal 
drop near 15 arcmin as well as the remaining EEB, EBB, and BBB 
signals below 5 arcmin can be explained with IA. This supports the¬ 
oretical models of IA that predict a insignificant contribution of IA 
correlations to the EEE signal above ~ 5 arcmin for CFHTLenS 
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dMerkel & S chafeil 1 201 4l : IValageasll2014l : ISemboloni et al.ll2008l) . 
In conclusion, further research is required to remove the remaining 
EEB, EBB, BBB systematics and to decide whether the tension be¬ 
tween data and theoretical models persist. Finally, F14 do not report 
a significant B-mode on scales below 5 arcmin for the second-order 
aperture statistics. But they agree with our finding of a EEB, EBB, 
and BBB signal on these scales for the third-order statistics. This 
suggests that the here reported PSF systematics become only rele¬ 
vant for higher-order cosmic shear statistics. 

As additional online material we provide a Monte-Carlo sam¬ 
ple of (fi m , erg) based on the posterior in the top left panel of Fig. 
flTland a set of 200 realisations of (Mf p ) data vectors that we pro¬ 
duced from our noisy clone simulations. 
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APPENDIX A: REFINED TEST FOR SHEAR 
SYSTEMATICS 

The level of residual systematics in the CFHTLenS lensing cata¬ 
logue is quantified by H12 using the two-point shear statistics. It is 
therefore not guaranteed that residual systematics for higher-order 
statistics would lead to the same field selection. For this paper, we 
therefore assess the level of systematics affecting three-point shear 
statistics by using a refinement of the H12 methodology. We apply 
the refined test to the 129 CFHTLenS pointings that have already 
passed the systematics criteria of HI2. Out of these 129 we reject 
another nine pointings for our cosmological analysis of the three- 
point shear statistics. The details of the refinement and its results 
for CFHTLenS follow below. 


Al Description of the method 

We adopt the residual systematics PSF model of H12: any resid¬ 
ual systematics is a linear combination of the stellar anisotropy e* a 
of all exposures a. Therefore, the observed galaxy ellipticity, e obs , 
originating from the lack of a perfect PSF correction, is given by: 

^exp 

e° bs = «i nt + 7< + Vi + , *=1,2. (Al) 

a=1 

This thus expresses, for each ellipticity component i, the observed 
(PSF-corrected) ellipticity e° ba as sum of (i) the intrinsic elliptic¬ 
ity ef“, (ii) the cosmic shear 7 ;, and (iii) the random noise r/i. In 
addition, (iv) the last term expresses the PSF residual systematic 
error at the location of the galaxy as a linear combination of the 
original PSF ellipticity e* directly measured from the stars in each 
exposure a out of n exp exposures; the a a are the coefficients of the 
linear combination. 

Moreover, following the conclusion from H12 we assume that 
the zero-lag correlations of the ellipticities already contain all the 
relevant information about the residual PSF correlations for the 
three-point correlation functions. For the PSF residuals described 
by the average measured zero-lag on a given pointing is 

/ obs obs obs\ / \ 1 / a / obs obs obs\\ / a 

(G <7 G, J = (7<7j7fc) + ej Gc )J , (A2) 
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where i,j,k = 1,2 indicate the projections of e obs and 'j, and 

(A(e° bs e° bs ef a )) (A3) 

n exp 

^ ^ OLaOibCX-c (ta,i€b,j tcjk') 


a,b,c= 1 


'‘exp 

+ 3 [ {jj'Yk) + (VjVk) + («r^ nt ) ] ° a ( e «.i) • 


The indices (a, b, c) indicate the various exposures, while the av¬ 
erage (...) is over all zero-lag triplets. Similarly, the PSF-galaxy 
cross-correlations are: 


/ obs obs ★ \ 

\ £i e i £a ’ k / 

= [ (7i7 j) + (ViVj) + (e mt e mt ) ] (C,k) 

n exp 

+ ^ ^ OLaOt-bOtc ^a,k^-b,i^-c,j ^ 

b,c= 1 


(A4) 


and 


n exp 

€a,j€b,k/ ^ ^ Q! c ^€ Cj ^b,k ^ • (A5) 

c =1 


We estimate the zero lag triplets by interpolating the stellar 
anisotropy at the position of the source galaxy. For the derivation of 
the equations {43HA5}, we assume that the PSF is uncorrelated 
with the intrinsic ellipticity, the shear and the random noise. In the 
absence of any systematics, equations (lA3T> (fA5ll have to vanish, 
hence a non-zero signal for either average can be used as indicator 
for systematics. However, in the presence of systematics the ex¬ 
pectation values of both dA3b and dA4b do directly depend on cos¬ 
mology through the terms ( 7 i 7 j 7 fc) and ( 7 i 7 y), and they depend 
on the details of IA through Therefore, for an evalua¬ 

tion of the significance of a non-zero signal we have to assume a 
fiducial model for the shear and IA correlations. Conversely, the 
expectation value of iA5l is free of these assumption so that we fo¬ 
cus on as systematics indicator in the following. Nevertheless 
some cosmology dependence enters for this indicator too because 
the variance of this indicator 

/ / obs * * \2\ 

£a,j£b,k) / = 

+ 

contains second-order correlations between 7 i and e^ nt on the right 
hand side; the intrinsic variance for the null hypothesis, i.e., no sys¬ 
tematics, has to be known to test for a presence of systematics. 
Contrary to d A 3 b and dA4b however, we expect a weak impact of 
the cosmology on the test results, though, as on small scales (zero 
lag) the variance ((e; nt + 7 ; + r/i) 2 ) ft; (( 4 nt ) 2 ) is likely domi¬ 
nated by shape noise. 


((<7* +7i +Vi) 2 ) ((C,j4,k) 2 ) 


“exp 

\ ^ * * * 

/ j &c€c,i€a,j€b,k 

c =1 


(A6) 


A2 Null hypothesis 

In the following, we describe the set of simulations we use 
for the measurement of the PDF of the systematics indicator 
(e° bs ejj£a,fc) f° r the null hypothesis. This distribution is the key 
element in the evaluation of the residual systematics for the 3- 
points function measurement. Using the same set of simulations 
described in H12, there are 160 realisations per pointing avail¬ 
able, each with a shear signal assigned from a projected mass map 


W3m2m1_y 



Figure Al. Top panel, comparison of the systematics indicators E i of 
W3m2ml_y (vertical black line) to the simulated distribution of the null 
hypothesis (histogram). The red solid line is the best-fitting Gaussian to 
the histogram. The Hi is consistent with no systematics. Bottom panel, the 
same as the top panel but for S 2 . Here the measured value is inconsistent 
with the null hypothesis on a high confidence level; the pointing is hence 
rejected for the cosmological analysis of 


extracted from independent lines-of-sights. These maps originate 
from the same N-body simulations that have been used for the 
clone. This means for the null hypothesis: 

• a WMAP5 cosmology; 

• no intrinsic alignments of the sources; 

• no correlations of intrinsic ellipticities and shear; 

• no B-modes of the shear field; 

• distribution of intrinsic shapes and source positions as in 
CFHTLenS. 

A shear value is assigned to each galaxy, depending on its 3D po¬ 
sition in the sky. Shape noise is added to the shear as an ellipticity 
component obtained by randomising the orientation of the elliptic¬ 
ity of the CFHTLenS galaxies. This procedure guarantees that the 
simulated catalogues have similar shape noise and intrinsic ellip¬ 
ticity characteristics as the real CFHTLenS pointing. Note that the 
simulated shape noise is marginally larger than in the real data be¬ 
cause the shear is also randomised when computing the noise, it is 
not removed; this is only a one percent effect and hence negligi¬ 
ble for the test. For each pointing, we create 1600 realisations by 
generating for each of the 160 lines-of-sight ten different noise re¬ 
alisations. The PSF ellipticit y at each galaxy lo cation, e* a is derived 
from the lensfA PSF model I lMiller et al.ll2013l) . For each pointing, 
the PDF of (e° bs £a,j £ a,fc) i s then constructed from the distribution 
of dA5b in the 1600 realisations. 

The probability of measuring a given value of 

n exp 

= > * = 1.2 (A7) 

a =1 

for a CFHTLenS pointing given its specific PSF and noise prop¬ 
erties can now be quantified against the assumptions of the null 
hypothesis. We obtain the zero lag correlation (e° bs £a,j e 6,fc) by in¬ 
terpolating e*j and e* k at the position of the galaxy with e obs . For 
simplicity, we restrict the analysis to the correlations (e° ha e* a ^e* a ^ 
measured in the same exposure a, and then sum over all the ex¬ 
posures of a given pointing; cross-correlations between exposures 
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Si S 2 null (Gaussian) 

|3i| < <n 90 (70%) 69 (53%) - 80(68.2%) 

|3ij < 2 ai 126(98%) 109(84%) - 123(95.4%) 

|Sij < 3<ri 129(100%) 124(96%) - 129(99.6%) 

Table Al. Number of pointings with |Sj| below <Ji (first line), 2ai (second 
line), and 3er.j (third line); ai is the standard deviation of the null hypoth¬ 
esis (no systematics). We indicate in parentheses the fractional value cor¬ 
responding to this number. We show the results both for Si (first column) 
and S 2 (second column) as well as the expectation for a null signal (third 
column). The total number of 129 pointings used here complies with the 
systematics criteria of HI2. 


are ignored. This strategy provides the strongest signal. We then 
judge the systematics significance of a value of Si in CFHTLenS 
by the correlation excess with respect to our simulated PDF. See 
Fig. IA1I based on the pointing W3m2ml_y, as an example. Here 
we plot the null hypothesis as histogram (black solid lines) that we 
successfully fit by a Gaussian distribution (red solid lines). The ver¬ 
tical solid lines are the corresponding actual values in CFHTLenS 
for comparison. We reject this pointing due to its large excess of 
S 2 . 


A3 Application to the CFHTLenS data 

We evaluate the cross-correlation defined in dA5b for each pair of 
exposures (a, b) and for different projections (i,j- k) of the vec¬ 
tors (e*, e obs ). From this, we estimate the statistics E; for every 
CFHTLenS pointing and test the values against the null hypothe¬ 
sis. For this test, we assume that a Gaussian distribution is a good 
approximation of the null distribution of Hi. This is a valid as¬ 
sumption for the two-points cross-correlations (see H12), and we 
have checked that it is also a valid assumption for histograms of 
Hi, which we fit by a Gaussian of mean and variance < 7 ; for ev¬ 
ery pointing. As expected, the averages iz, are always consistent 
with zero. 

As a null E; is supposed to obey Gaussian statistics, 31.80% 
of the pointings should have |E;| > cr;, 4.6% |E;| > 2cr;, and 
0.04% | Si | > 3(7;. The Table I A ll compares systematics indicators 
of 129 CFHTLenS pointings to the null hypothesis with respect to 
the < 7 ;, 2cr i, and 3(7; thresholds. These 129 pointings have been pre¬ 
selected by the criteria defined in HI2. Within 3(7; the statistics of 
| Ei | is consistent with the null hypothesis, whereas | S 2 1 reveals too 
many outliers with | S 2 1 > 3ct;. 

For the final cosmological analysis, we decided to re¬ 
ject pointings that are within the Gaussian 1% tail of the 
null hypothesis. The false-positive rate of our test is con¬ 
sequently 1%. Based on this cut, we reject the following 
pointings: Wlm3p3_i, Wlp3pl_y, W2mlm0_i, W3m3m0_i, 
W3m2ml_y, W3plml_i, W4m3pl_i, W4m3m0_i, and 
W4m3p3_y. All these fields are rejected owing to too high values 
of | S 2 1 alone. The signal with and without the rejected fields can 
be seen in the left panel of Fig. [4] 
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